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ABSTRACT 


A  combined  experimental  and  analytical  program  has  been  conducted  to 
examine  transitional,  accelerating  boundary  layer  flows  with  high  levels  of 
freestream  turbulence.  This  program  was  designed  to  complement  a  parallel 
study  conducted  previously  under  AFOSR  funding  (Contract  No.  F49620-78-C- 
0064).  The  earlier  program  focused  on  measurement  of  transitional  heat 
transfer  distributions  for  four  combinations  of  streamwise  acceleration  and 
freestream  turbulence.  The  present  program  was  designed  to  document  the 
boundary  layer  turbulence  structure  and  spectral  distributions  for  the  same 
four  test  conditions.  The  results  from  the  present  program  have  shown  that 
transition  in  accelerating  flows  consists  of  an  acceleration  dominated  stage 
of  slowly  developing  intermittency  followed  by  a  second  stage  with  the  same 
general  characteristics  as  zero-pressure-gradient  transition.  Conditionally 
sampled  fluctuating  velocity  profile  measurements  indicated  that  the  boundary 
layer  turbulence  was  highly  anisotropic  in  the  early  stages  of  transition. 
Conditionally  sampled  mean  velocity  measurements  showed  that  within  the 
intermittent  turbulent  patches  the  mean  velocity  profiles  were  very  similar  to 
those  of  an  equilibrium  turbulent  boundary  layer.  Spectral  distribution  data 
indicated  that  preferred  amplification  of  the  most  unstable  (as  predicted  by 
linear  stability  theory)  frequencies  occurred  upstream  of  the  onset  of 
transitional  bursting. 

In  addition  to  the  experimental  portion  of  this  investigation,  numerical 
experiments  were  undertaken  to  assess  the  ability  of  currently  existing 
methods  to  predict  heat  transfer  during  transition  in  accelerating  flows. 
Comparisons  of  these  numerical  results  with  the  present  experimental  data 
indicate  that  excellent  predictions  of  both  heat  transfer  and  boundary  layer 
development  can  be  achieved  with  case-specific  knowledge  of  the  streamwise 
intermittency  distribution.  General  predictions,  which  required  intermittency 
correlations  as  input,  were  much  less  satisfactory. 
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I.  INTRODUCTION 


Accurate  calculation  of  airfoil  boundary  layer  development  is  a  prerequi¬ 
site  for  the  prediction  of  turbine  heat  transfer  distributions  and  aerodynamic 
performance  (efficiency).  Such  calculations,  however,  are  exceptionally  dif¬ 
ficult  in  that  the  nature  of  boundary  layers  on  turbine  airfoils  and  the 
numerous  mechanisms  that  affect  them  are  extremely  complicated.  Turbine  air¬ 
foil  boundary  layers  have  regions  where  they  are  laminar,  transitional 
(forward  and  reverse)  and  fully  turbulent.  They  are  also  subject  to  the 
interacting  effects  of  streamwise  curvature,  freestream  turbulence,  pressure 
gradients  and  three-dimensionality. 

Precise  prediction  of  airfoil  profile  heat  loads  is  of  critical  impor¬ 
tance  because  the  maximum  turbine  inlet  temperatures  nearly  always  exist  in 
the  raidspan  region  of  the  airfoil  away  from  the  endwalls.  Fortunately,  for 
the  vast  majority  of  turbine  configurations,  the  airfoil  profile  boui.dary 
layers  are  two-dimensional  for  most  of  the  span.  Only  near  the  root  and  tip 
is  this  situation  complicated  by  the  presence  of  three-dimensional  secondary 
endwall  flows.  There  is,  then,  an  important  need  to  be  able  to  accurately 
compute  two-dimensional  boundary  layer  development  and  heat  transfer  for  the 
turbine  environment. 

The  boundary  layers  which  develop  along  the  suction  and  pressure  surfaces 
of  modern  high  camber  turbine  airfoils  have  vastly  different  streamwise  pres¬ 
sure  histories.  Along  the  suction  (convex)  surface  there  is  typically  a  half 
chord  of  moderately  high  acceleration  upstream  of  the  airfoil  minimum  pressure 
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region  followed  by  a  moderate  deceleration  to  the  airfoil  trailing  edge. 
Depending  upon  the  fore-chord  acceleration  level  and  the  freestream  turbulence 
the  suction  surface  boundary  layer  may  become  transitional  anywhere  from  up¬ 
stream  of  the  maximum  velocity  region  to  as  far  downstream  as  the  adverse 
pressure  gradient  zone. 

Along  the  pressure  (concave)  surface  the  developing  boundary  layer  exper¬ 
iences  approximately  half  a  chord  of  low,  near-constant  velocity  followed  by  a 
length  of  extremely  high  acceleration  to  the  airfoil  trailing  edge.  The  state 
of  turbulence  in  the  pressure  surface  boundary  layer  is  determined  by  an 
extremely  complex  interaction  of  phenomena.  In  addition  to  the  destabilizing 
effects  of  the  high  freestream  turbulence,  the  concave  wall  curvature  produces 
destabilizing  streamwise  vortex  systems  within  the  boundary  layer.  Countering 
these  turbulence  promoting  effects,  the  extreme  aft-chord  accelerations  can  be 
sufficient  to  relaminarize  even  fully  turbulent  boundary  layers.  Since  local 
heat  transfer  rates  are  very  sensitive  to  the  state  of  the  boundary  laver 
(laminar,  transitional  or  turbulent)  the  ability  to  predict  transition  onset 
and  length  is  very  important.  As  an  illustration  of  the  importance,  it  is  not 
uncommon  for  the  boundary  layer  to  be  transitional  for  over  half  the  chord  for 
a  typical  turbine  airfoil. 

There  is  a  widely  recognized  need  to  improve  existing  analytical  models 
for  prediction  of  boundary  layer  development  (profile  loss)  and  heat  transfer 
distributions  for  transitional  turbine  airfoils.  One  family  of  cases  for 
which  particularly  large  discrepancies  between  prediction  and  experiment  have 
been  observed  are  transitional  flows  with  streamwise  acceleration  and  high 
freestream  turbulence.  Brown  and  Martin  (Ref.  1)  made  special  note  of  the 
current  lack  of  adequate  understanding  of  this  problem  in  their  comprehensive 


heat  transfer  review  paper.  Additional  examples  of  important  discrepancies 
between  state-of-the-art  airfoil  boundary  layer  predictive  schemes  and  the 
best  available  experimental  data  are  given  by  Han,  Chait,  Boyee  and  Rapp  (Ref. 
2),  Hylton,  Mihelc,  Turner,  Nealy  and  York  (Ref.  3)  and  Rae,  Taulbee, 

Civenskas  and  Dunn  (Ref.  4). 

Currently  the  most  highly  mature  portion  of  the  analytical  prediction  of 
transition  is  the  computation  of  unstable  normal  modes  in  boundary  layers  with 
small  (linear)  single  frequency  (Tollmien-Schl icht ing)  disturbances.  Much 
effort  has  also  been  expended  to  develop  methods  for  predicting  the  amplifica¬ 
tion  process  through  which  the  single  frequency  disturbances  grow  to  nonlinear 
levels  and  culminate  in  turbulent  spot  formation  (see  Ref.  5  for  a  description 
of  the  e9  and  "modified  e9"  amplitude  ratio  methods).  In  addition,  Mack  (Ref. 
6)  has  developed  techniques  for  applying  stability  theory  to  boundary  layers 
exposed  to  a  wideband  disturbance  spectrum  (e.g.,  freestream  isotropic  or 
near- isotropic  turbulence).  Using  numerical  computational  techniques  for  a 
family  of  Falkner-Skan  family  acceleration  profiles,  Mack  has  successfully 
predicted  the  trends,  known  long  from  experimental  testing,  of  streamwise 
acceleration  to  retard  and  of  freestream  turbulence  to  promote  the  transition 
process  . 

The  analytical  techniques  discussed  above  apply  exclusively  to  low- 
intensity  disturbances  in  laminar  boundary  layers.  For  large  levels  of  exter¬ 
ior  disturbances  both  Morkovin  (Ref.  5)  and  Reshotko  (Ref.  7)  subscribe  to  the 
concept  of  the  "high  intensity  bypass"  (complete  jump  of  the  linear  amplifica¬ 
tion  stage  of  transition).  The  bypass  mode  of  transition  is  thought  to  occur 
when  disturbances  to  the  boundary  layer  are  so  large  that  they  are  already  at 
non-linear  levels  without  Tol lmien-Schl icht ing  amplification.  The  available 


experimental  evidence  supports  this  concept  in  that  the  highest  freestream 
turbulence  levels  for  which  Tollmien-Schlicht ing  waves  have  ever  been  detected 
are  0.45%  for  flat  wall  transition  (Ref.  8)  and  0.5%  for  turbine  airfoil  tran¬ 
sition  (Ref.  9  and  10).  One  portion  of  this  present  program  was  designed  to 
provide  new  information  on  the  "bypass"  mode.  For  relatively  high  freestream 
turbulence  levels  the  spectral  distribution  of  the  boundary  layer  disturbances 
upstream  of  the  onset  of  transition  were  examined  to  detect  preferential  amp¬ 
lification  of  the  bandwidth  of  predicted  unstable  frequencies. 

Downstream  of  the  "onset"  of  transition,  i.e.,  downstream  of  the  first 
appearance  of  turbulent  bursts  in  the  test  boundary  layers,  the  present  study 
examined  the  development  of  the  intermittent  boundary  layer  turbulence  struc¬ 
ture.  A  common  method  of  dealing  with  transition  in  current  finite-difference 
boundary  layer  computations  is  to  presume  that  it  occurs  over  zero  length. 

The  flow  is  assumed  to  be  fully  laminar  at  one  computation  station  and  to  have 
an  equilibrium  turbulence  structure  at  the  next.  Other  analytical  models  have 
attempted  to  incorporate  the  experimentally  observed  facts  that  (1)  real  tran¬ 
sition  extends  over  a  distance  of  many  boundary  layer  thicknesses  (sometimes 
hundreds)  and  (2)  the  flow  in  the  transitional  region  is  intermittent  in 
character — sometimes  laminar  and  sometimes  turbulent.  Some  of  the  more  widely 
used  transitional  flow  turbulence  models  of  this  type  include  those  of 
McDonald  and  Fish  (Ref.  10),  Price  and  Harris  (Ref.  11  and  12)  and  Forest 
(13).  All  of  these  models  rely  upon  the  assumption  that  at  any  location  along 
the  length  of  transition  certain  turbulence  quantities  (e.g.,  the  Reynolds 
stress  structural  coefficient  u'v'/q^  for  Ref.  10  or  the  eddy  viscosity  for 
Ref.  13)  are  equal  to  the  equilibrium  turbulent  boundary  layer  value  times  the 


local  interraittency .  Since  all  these  models  rely  upon  unverified  assumed 
distributions  of  turbulence  quantities  they  can  be  accurately  described  as 
speculative  in  nature.  The  program  was  designed  to  provide  the  data  required 
to  incorporate  more  realistic  physics  into  such  transition  models. 

This  present  program  was  designed  to  complement  a  parallel  study 
conducted  previously  under  AFOSR  funding  (Contract  No.  F49620-78-C-0064 ) .  The 
earlier  program  focused  on  measurement  of  transitional  heat  transfer 
distributions  for  four  combinations  of  streamwise  acceleration  and  freestream 
turbulence.  The  present  program  was  designed  to  document  the  boundary  laver 
turbulence  structure  for  the  same  four  test  conditions.  The  measurements  from 
’"'ife  earlier  program  combined  with  the  present  results  provide  a  comprehensive 
examination  of  accelerating  transitional  boundary  layer  flow.  Analytical 
results  include  comparisons  between  these  various  measurements  and  a  variety 
of  numerical  transitional  boundary  layer  computation  procedures. 


II.  DESCRIPTION  OF  TEST  EQUIPMENT 


1.  UTRC  BOUNDARY  LAYER  WIND  TUNNEL 

All  experimental  data  for  the  present  investigation  were  obtained  in  the 
United  Technologies  Research  Center  (UTRC)  Boundary  Layer  Wind  Tunnel.  A 
complete  description  of  this  facility  is  given  in  Reference  14.  This  tunnel 
was  designed  for  conducting  fundamental  studies  of  two-dimensional,  incompres¬ 
sible  flat  wall  boundary  layer  flow.  Incorporated  in  the  tunnel  is  a  versa¬ 
tile,  adjustable  test  section  constructed  so  that  laminar,  transitional,  or 
turbulent  boundary  layers  can  be  subjected  to  favorable,  zero,  or  adverse 
pressure  gradients.  In  addition,  test  boundary  layers  can  be  subjected  to  a 
wide  range  of  freestreara  turbulence  levels.  Low  freestream  turbulence  flows 
can  be  investigated  in  this  facility  since  it  is  designed  to  have  a  very  low 
residual  test  section  turbulence  level.  Higher  turbulence  levels  can  be  gen¬ 
erated  within  the  test  section  through  the  use  of  various  rectangular 
grids. 

An  overall  sketch  of  the  Low  Speed  Boundary  Layer  Tunnel  is  shown  in 
Figure  1.  The  tunnel  is  of  recirculating  design  and  consists  of  a  blower,  a 
settling  chamber /p lenum,  a  contraction  nozzle,  the  boundary  layer  test  sec¬ 
tion,  a  downstream  diffuser,  and  a  return  duct.  The  settling  chamber /plenum 
consists  of  a  series  of  perforated  part  span  baffles  which  even  out  gross 
irregularities  in  the  flow  from  the  blower  and  a  honeycomb  which  removes 
large-scale  flow  swirl.  Downstream  of  the  honeycomb  are  a  series  of  fine  mesh 
damping  screens  which  progressivly  reduce  both  the  flow  nonuniformity  and  the 
residual  tunnel  turbulence  level.  A  nozzle  with  a  2.8:1  contraction  ratio 


mounted  downstream  of  the  damping  screens  accelerates  the  flow  to  produce  the 
required  test  section  Reynolds  number.  Following  the  contraction  nozzle  the 
flow  passes  through  the  34-in.  wide  flat  wall  boundary  layer  test  section.  At 
the  entrance  to  the  test  section  an  upstream  facing  scoop  bleed  assembly  forms 
the  leading  edge  of  the  boundary  layer  test  surface.  The  purpose  of  this 
leading  edge  bleed  scoop  is  to  remove  all  the  flow  near  the  tunnel  upper  wall. 
With  this  arrangement  the  test  section  flow  consists  of  the  uniform  "core" 
flow  from  the  main  contraction  nozzle.  The  scoop  assembly  consists  of  a  two- 
stage  leading  edge  adjustable  bleed  attached  to  the  flat  wall  boundary  layer 
test  surface.  The  upstream  and  by  far  the  larger  of  the  two  scoops  removes 
the  flow  nearest  the  upper  wall  of  the  contraction  exit  duct.  This  large 
scoop  is  intended  to  trap  both  the  two-dimensional  boundary  layer  which 
develops  along  the  contraction  nozzle  wall  and  the  vortices  which  develop  in 
the  contraction  corners.  The  local  scoop  flow  rate  can  be  adjusted  to  produce 
uniform  pressure  (in  the  transverse  direction)  at  the  static  taps  along  the 
entire  scoop.  The  downstream  and  much  smaller  of  the  two  scoops  is  mounted 
directly  on  the  front  edge  of  the  boundary  layer  test  surface.  The  test  sec¬ 
tion  boundary  layer  begins  growing  at  the  leading  edge  of  this  smaller  scoop. 
The  leading  edge  of  the  small  downstream  scoop  is  a  4  x  1  elliptical  shape  in 
order  to  prevent  a  local  separation  bubble  and  a  premature  transition  of  the 
test  surface  boundary  layer. 

The  test  section  of  the  Boundary  Layer  Wind  Tunnel  consists  of  a  flat, 
smooth,  aluminum  upper  wall  test  surface,  a  lower  flexible,  adjustable  stain¬ 
less  steel  wall  and  transparent  vertical  sidewalls.  The  vertical  sidewalls 
were  constructed  of  plexiglass  to  facilitate  positioning  of  boundary  layer 
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probes  and  for  purposes  of  conducting  flow  visualization  studies.  Downstream 
of  the  test  section  a  d i f fuser/corner  combination  reduces  the  flow  velocity 
and  delivers  the  flow  to  the  return  duct.  Mounted  in  this  return  duct  are  an 
air  filter  and  a  liquid  chilled  heat  exchanger  which  controls  and  stabilizes 
the  tunnel  air  temperature  at  approximately  70°F. 
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2.  TURBULENCE  GENERATING  GRIDS 

As  described  in  Reference  14,  this  wind  tunnel  has  a  relatively  low  resi¬ 
dual  test  section  turbulence  level  (<  1/4%).  Higher  turbulence  levels  re¬ 
quired  for  this  study  were  generated  by  inserting  various  square  array  biplane 
grids  constructed  from  rectangular  bars  at  the  entrance  to  the  main  tunnel 
contraction  (see  Fig.  1).  Three  turbulence  generating  grids  were  designed 
using  the  correlations  of  Reference  15.  The  grids  will  be  referred  to  as 
Grids  1,  2,  and  3  corresponding  to  mesh  widths,  M,  of  7/8,  2  9/16,  and  7  in. 
(see  Fig.  2).  Details  of  the  grid  configuration  are  given  in  Reference  14. 
This  present  arrangement  differs  from  that  used  for  nearly  all  the  earlier 
investigations  of  this  subject  in  which  the  turbulence  grids  were  located  in 
the  test  section  just  upstream  of  the  boundary  layer  test  surface.  The  bene¬ 
fits  derived  from  locating  the  grids  at  the  contraction  entrance  were  that  the 
generated  turbulence  was  more  homogeneous  and  had  a  lower  decay  rate  along  the 
test  section.  Since  grid  generated  turbulence  decays  approximately  as  u'/U  a 
(x/b)-5/7  (Ref.  15),  the  change  in  turbulence  level  with  distance  along  the 
test  section  was  reduced  by  increasing  the  distance  from  the  grid  to  the  test 
section  entrance.  In  addition,  the  results  of  Reference  15  indicate  that 
approximately  10  grid  mesh  lengths  are  required  to  establish  a  uniform 
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turbulent  flow.  Locating  the  grid  a  distance  upstream  of  the  test  section 
requires,  of  course,  a  more  coarse  grid  to  achieve  a  given  test  section  tur¬ 
bulence  intensity. 

Another  effect  considered  was  the  expected  influence  of  the  contraction 
on  the  components  of  the  grid  generated  turbulence.  It  was  recognized  that 
rearrangement  of  the  relative  magnitudes  of  the  turbulence  components  would 
occur  due  to  the  contraction.  However,  since  the  contraction  ratio  was  small 
(2.8),  it  was  concluded  that  any  effects  of  induced  anisotropy  would  be  small 
in  comparison  to  the  advantages  gained  in  homogeneity  and  reduced  decay  rate. 
As  part  of  an  earlier  study  of  transitional  heat  transfer  in  this  same  facil¬ 
ity  (Refs.  16  and  17)  all  three  components  of  the  test  section  tubulence  were 
documented  for  these  same  test  grids.  Repeat  measurement  (at  a  reduced  number 
of  locations)  of  these  three-component  freestreara  turbulence  data  were  also 
obtained  as  part  of  this  present  program. 

3.  TEST  SECTION  INSERTS 

Two  test  section  inserts,  to  be  installed  opposite  the  flat  aluminum  test 
wall,  were  designed  to  produce  flows  with  a  nearly-constant  acceleration  par¬ 
ameter  (K  =  v/u2  3U/3X  =  0.2  and  0.75  x  10-6,  wedges  1  and  2  respectively). 
These  two  test  section  inserts,  sketches  of  which  are  shown  in  Figure  3,  con¬ 
sist  of  simple  wedge-shaped  bodies  with  a  2  in.  long  steeper  wedge  attached  to 
the  leading  edge.  These  "modified  shape"  wedges  were  designed,  using  the 
inviscid  potential  flow  analysis  of  Reference  18  to  provide  a  near  constant 
acceleration  of  the  test  section  flow  along  the  entire  test  wall.  Important 
dimensions  of  the  two  wedge  inserts  are  given  in  Figure  3.  Probe  traverse 
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slots,  aligned  with  slots  in  the  test  section  bottom  wall,  have  been  incorpor¬ 


ated  into  the  wedge  inserts. 


4.  INSTRUMENTATION  AND  DATA  ACQUISITION 


a.  Hot  Wire  probes — Measurements  of  fluctuating  velocities  in  the  test 


boundary  layers  were  obtained  using  multi-element  hot  wire  aneraometry 


techniques.  Both  vertical  and  horizontal  x-type  two  wire  probes  as  well  as 


single  horizontal  wire  probes  were  employed.  In  order  to  minimize  potential 


errors  for  these  measurements  (errors  largely  arising  from  the  inherent  mean 


velocity  gradients  in  the  flows  and  the  finite  probe  size)  the  hot  wire  probes 


were  custom  designed  and  fabricated  specifically  for  boundary  layer 


measurements.  A  comprehensive  description  of  the  design  of  these  probes  is 


given  in  Reference  19.  Included  in  Reference  19  are  detailed  examples  of 


velocity  statistics  measured  in  zero  pressure  gradient  boundary  layers. 


Following  a  preliminary  "burn  in"  and  a  resistance-temperature  calibra¬ 


tion  each  sensor  was  calibrated  for  velocity  and  angular  sensitivity  in  a  low- 


turbulence  1  1/2-in.  dia.  jet  flow  for  approximately  twenty  jet  flow  speeds 


ranging  from  7  to  130  ft/s.  The  mean  response  equation  of  each  sensor  was 


assumed  to  be  of  the  form 


Nu  =  A1  +  BlRe 


which  can  be  algebraically  manipulated  to 


A2(Rs  ♦  R„)2  „  .  .  B2(Rs  *  «J2 


T0-76  (Tw  -  T)  +  - - —  (Tw  -  T)  Ue0-45 

Rw 


where  Ew  =  wire  voltage 


R  =  probe  body,  cable  and  internal  anemometer  resistance 


R.,  =  sensor  resistance 
w 

T  =  air  temperature 
Tw  =  sensor  temperature 
UE  =  effective  velocity 
^2’  ®2  =  etnpi-rical  constants 

The  constants  A2  and  B2  were  determined  for  each  sensor  from  a  least- 
squares  fit  of  the  data  to  Equation  (2).  Next,  for  the  x-wire  probes,  pitch 
angle  versus  voltage  data  were  obtained  with  the  probes  rotated  from  +20°  to 
-20°  in  steps  of  5°.  The  center  of  pitch  coincided  with  the  intersection  of 
the  wires  of  the  x.  The  angular  sensitivity  of  the  wires  was  assumed  to 
conform  to  Champagne's  k2  law  (Ref.  20), 

Ug 2 ( 0 )  =  U2 (0  =  0)  (cos2  0  +  k2 s in20 )  (3) 

where  0  =  angle  between  wire  and  direction  normal  to  the  flow  (0  *  45° 
with  probe  stem  normal  to  the  flow) 

Ug  =  effective  velocity 

Using  a  least-squares  routine  to  find  a  best  fit  of  the  pitch-voltage 
data  to  Champagne's  equation,  optimum  values  of  k  were  determined  for  each 
sensor. 

In  summary,  the  temperature-resistance,  mean  velocity  and  pitch  calibra¬ 
tions  were  used  to  determine  the  following  calibration  constants. 

(1)  R^2  “  sensor  resistance  at  32°F 

(2)  a  -  temperature-resistance  coefficient 

(3)  A2  and  B2  -  empirical  constants  (Eq.  2) 

(4)  k  -  empirical  constant  (Eq.  3). 

b,  Anemometry  and  signal  processing  equipment--Thermo  Systems  Inc. 

(TSI)  model  1050  constant  temperature  anemometers  were  employed  to  drive  the 
hot  wire  sensors.  Standard  square-wave  techniques  were  used  to  assure  that 
the  sensor /anemometer  frequency  response  exceeded  50  KHz.  Signals  from  the 
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anemometers  were  first  passed  through  5  KHz  low-pass  anti-aliasing  filters, 
reduced  with  a  precise  D.  C.  offset,  passed  through  a  wide  band  amplifier 
(Preston  Model  8300  XWB)  and  then  digitized  using  a  TSI  Model  1075  Multichan¬ 
nel  Digitizer.  A  feature  of  this  particular  analog-to-digitial  converter 
which  is  important  to  this  application  is  that  the  various  channels  are  sam¬ 
pled  and  held  simultaneously.  This  simultaneous  sample-hold  feature  permits 
cross-product s  of  the  various  fluctuating  quantities  to  be  computed.  The  buck 
and  gain  system  permits  the  maximum  12  bit  resolution  of  the  digitizer  to  be 
util ized . 

Signals  were  sampled  for  5  seconds  at  7813  samples/sec  for  a  total  of 
39060  continuous  numbers  per  channel.  The  bandwidth  of  frequency  resolution, 
then,  was  1/5  to  3900  Hz.  A  Perkin  Elmer  3205  minicomputer  with  800  K  Bytes 
available  RAM  was  used  to  control  the  digitizer  and  temporarily  store  a  stream 
of  data.  Longer  term  storage  of  the  raw  voltage-time  records  was  achieved 
using  25  M  Byte  hard  disks  and  digital  tape.  Reduction  of  the  voltage-time 
records  to  velocity-time  and  mean  velocity  statistics  was  also  accomplished 
using  the  PE  3205.  Computed  quantities  included  the  mean  values  of  first 
through  fourth  moments  of  the  velocity  fluctuations  and  the  double  (Reynolds 
stress)  and  triple  velocity  products. 
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III.  TEST  CONDITIONS 


1.  TEST  SECTION  VELOCITY  DISTRIBUTIONS 

The  acceleration  configurations  (1  and  2)  were  designed  to  produce  "sinl 
or  constant  acceleration  parameter  flows.  Idealized  "sink"  flow  is 
represented  in  the  following  sketch. 


^ZZZ/Z/ZZZ/ZZZZ^l 


where  X  is  the  distance  from  the  beginning  of  acceleration  and  Xg  is  the 
distance  to  the  potential  "origin".  The  velocity  distribution  for  idealized 
sink  flow  is 

u  =  c(Xs  -  50-1  (a 

The  velocity  distribution  produces  a  constant  acceleration  parameter 

K  *  v/U^  3U/3J?  *  v/C  ■  constant  (5 

For  the  present  test  program  the  wedges  were  carefully  adjusted  to 
reproduce  the  exact  velocity  distributions  employed  in  References  16  and  17. 
In  this  way  the  present  intermittent  turbulence  structure  data  and  the  wall 
Stanton  numbers  obtained  in  the  previous  study  can  be  combined  to  form  a 
single  comprehensive  set  of  transitional  boundary  layer  information.  The 


it 


velocity  distributions  reported  in  References  16  and  17  are  reproduced  in 


Figure  4  and  6.  The  present  data  agreed  within  -  1%  of  these  previously 


measured  distributions.  An  examination  of  Figures  4  and  5  reveals  that  the 


flows  generated  in  the  two  acceleration  configurations  closely  approximate 


ideal  sink  flow.  Test  boundary  layers  developing  along  the  flat  test  wall 


were  subjected  to  nearly  constant  acceleration.  The  velocity  ratio  curve 


plotted  for  Figure  4  is  the  analytical  relationship  determined  from  the 


velocity  distribution  data: 


U  =  14,800  (200  -  X)-1 -066 


where  U  is  in  ft/s  and  X  is  in  inches.  The  acceleration  parameter  K  was 


calculated  using  this  same  analytical  relationship.  As  shown  in  Figure  4  the 


acceleration  parameter  is  nearly  constant  along  the  entire  test  wall  with  a 


value  very  near  the  design  target  of  K  =  v/U2  3U/  <$X  =  0.2  x  10_o. 


For  Figure  5  the  velocity  ratio  and  acceleration  plots  were  calculated 


from  the  relationship 


U  =  3730  (83.3-X)'1 -075 


where  again  U  is  in  ft/s  and  X  is  in  inches.  The  streamwise  acceleration 


parameter  was  nearly  constant  along  the  entire  test  wall  with  K  =  0.76  x 


10_u  . 
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2.  FREESTREAM  TURBULENCE  DISTRIBUTIONS 


Boundary  layer  turbulence  statistics  were  measured  at  various  streamwise 
stations  along  the  test  wall  for  the  following  four  acceleration/turbulence 
comb inat ions : 


Flow 

Condition 

K 

x  10° 

Grid 

No. 

Denoted 

1 

0.20 

1 

K.20G1 

2 

0.20 

2 

K.20G2 

3 

0.75 

2 

K.75G2 

4 

0.75 

3 

K.75G3 

These  combinations  reproduced  the  same  test  conditions  for  which  mean 
velocity  and  temperature  profiles  and  wall  Stanton  number  were  recorded  for 
References  16  and  17. 

As  part  of  References  16  and  17  the  axial  distributions  of  the  freestream 
turbulence  intensity  and  streamwise  integral  length  scale  were  documented  for 
the  four  test  cases.  All  three  components  of  the  freestream  turbulence  were 
recorded  at  numerous  streamwise  stations.  These  previous  freestream 
turbulence  measurements  were  obtained  using  x-type  cylindrical  hot  film  probes 
and  linearized  constant  temperature  anemometers.  The  individual  velocity 
components  were  separated  using  sum-and-d i f f erence  circuits.  Figure  6 
presents  the  turbulence  intensity  and  streamwise  length  scale  data  obtained  at 
the  test  section  entrance  for  Reference  16.  These  data,  obtained  for  all 
three  turbulence  generating  grids,  indicate  that  both  the  intensity  and  scale 
were  independent  of  flow  speed  over  the  range  of  testing.  Multi-component 


freestream  turbulence  intensity  data  obtained  for  the  present  program  outside 
the  boundary  layers  are  presented  in  Figures  7  through  12.  These  new  data, 
obtained  with  the  previously  described  hot  wire  and  A/D  data  acquisition 
system,  are  plotted  along  with  the  earlier  data  of  Reference  16.  The  three 
components  of  turbulence  measured  for  K  =  0.2  x  10" 6  and  Grids  1  and  2  are 
presented  in  Figures  7  and  8  respectively.  The  present  results  are  in  very 
good  agreement  with  the  earlier  measurements.  Agreement  between  the  present 
and  previous  results  is  also  excellent  for  Figure  9  in  which  the  total 
turbulence  distributions  for  the  two  cases  are  compared.  The  individual 
component  distributions  for  Grids  2  and  3  with  K  =  0.75  x  10-°  (Figs.  10  and 
11)  compared  less  favorably  to  the  similar  previous  data.  In  contrast, 
however,  the  total  turbulence  measured  for  these  two  cases  (Fig.  12)  agreed 
very  well  with  the  earlier  results.  The  apparent  disagreement  shown  in  Figure 
10  and  11  resulted  from  the  fact  that  the  two  sets  of  intensity  measurements 
were  obtained  at  different  distances  from  the  wall.  The  data  of  Reference  16 
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IV.  CONDITIONAL  SAMPLING  TECHNIQUES 


The  detection  of  intermittent  turbulent  bursts  in  the  test  boundary  layer 
flows  was  accomplished  using  conditional  sampling  techniques.  These  proce¬ 
dures  are  well  established  as  an  accurate  means  for  discrimination  of  instan¬ 
taneous  flow  conditions  and  are  thoroughly  described  in  References  22  through 
24.  In  brief,  a  fluctuating  velocity  record  was  numerically  examined  to  de¬ 
termine  which  portions  should  be  designated  as  turbulent.  Specifically,  the 
conditional  sampling  criteria  were  selected  such  that  periods  of  high  fre¬ 
quency,  larger  magnitude  fluctuations  of  velocities  and/or  Reynolds  stresses 
were  identified  as  turbulent.  The  turbulent  criterion  were  met  when  a  detec¬ 
tor  function  (averaged  over  a  very  short  smoothing  period)  exceeded  a  pre¬ 
determined  threshold  level.  The  particular  choices  for  the  detector  function, 
smoothing  period  and  threshold  level,  then,  each  influenced  the  turbulent/non- 
turbulent  discrimination  process. 

Three  independent  detector  functions  were  examined  as  part  of  this  study. 

For  the  single  horizontal  hot  wire  records  the  square  of  the  streamwise  veloc- 
3u  2 

ity  fluctuations  — )  was  selected.  Squares  of  the  Reynolds  fluctuation 
3nC*  ("ft”)  were  used  for  the  vertical  and  horizontal  wires  respect¬ 
ively.  As  demonstrated  in  References  22  through  24,  use  of  the  square  of  the 
fluctuations  provides  a  clearer  discrimination  than  the  temporal  derivative 
alone . 
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Using  the  instantaneous  values  of  these  three  detector  functions 
velocity-time  and  Reynolds  stress-time  records  were  examined,  as  follows,  to 
empirically  select  smoothing  times  and  threshold  levels.  As  described  unh-r 
Data  Acquisition,  the  P-E  minicomputer  was  employed  to  reduce  the  x-tvpe  or 
single-wire  signals  to  velocity-time  records.  A  comprehensive  signal  proces¬ 
sing  software  package  (SPAG  -  Cranfield  Data  Systems)  then  permitted  conven¬ 
ient  displays  of  the  velocity-time  records  on  a  computer  monitor.  As  an 
example  the  u,  v  and  uv  product  over  a  period  of  time  could  be  displayed  on 
the  screen.  The  conditional  sampling  discrimination  procedure  was  programmed 
to  produce  a  display  showing  times  when  the  turbulent  criterion  was  satisfied. 
Hundereds  of  sample  velocity  records  covering  the  entire  test  ranges  of  free- 
stream  turbulence  and  intermit tency  were  then  tested  against  candidate  smooth¬ 
ing  times  and  threshold  levels.  The  values  eventually  selected  proved  to  be 
capable  of  reliably  identifying  turbulent  periods  over  a  very  broad  range  of 
conditions.  A  smoothing  period  of  0.0005  sec.  was  selected  for  use  with  all 
test  records.  This  smoothing  window  corresponds  to  approximately  60  -  100 
Kolmogorov  length  scales  for  most  of  the  data.  Very  similar  smoothing  periods 
were  employed  for  both  References  22  and  23.  The  threshold  level  selected 
consisted  of  a  multiple  of  an  upstream  reference  (pre-bursting,  near-wall) 
value  of  that  quantity  for  that  particular  test  boundary  layer.  For  the 

uv  and  uw  records  the  criteria  for  the  beginning  of  bursts  was  that  ( ~~~)  "> 

.  -  f3uv%2  3uw->2  „  ,  „  /Suw-,2  .  ,  _  ,  .  , 

10  I— — J  or  — — )  >  10  [— — J  respectively.  For  the  single 

v3t  'upstream  3t  '  v3t  'upstream  v  j  & 

wire  a  burst  beginning  was  satisfied  if  (-i-rl  >  6 

6  6  '•3t'  1 3t  -'ups  t  ream 

Identification  of  the  end  of  turbulent  periods  proved  to  be  considerably 
more  difficult.  The  leading  edges  of  bursts  were  consistently  sharp  (hence 
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the  large  multiplication  factor  on  the  threshhold  criterion).  There  were, 
however,  usually  periods  within  a  burst  where  the  fluctuations  were  relatively 
reduced  but  still  very  large  compared  to  the  upstream  "reference"  value. 
Examination  of  hundreds  of  velocity-time  displays  revealed  that  the  end  of 
bursts  consisted  of  a  "tailing  off"  of  the  high-frequency  fluctuations  by  a 
return  to  "upstream  reference"  conditions.  These  observations  were  accommo¬ 
dated  into  the  conditional  sampling  routine  by  requiring  that  8  smoothing 
periods  in  a  row  fall  below  a  "turbulent-end"  criterion  before  a  "stop"  was 
identified.  The  possible  error  introduced  by  this  logic  was  that  2  or  more 
short  bursts  could  be  lumped  together  into  a  single  artificially  long  burst. 
Ensemble  averaged  velocity  profile  data  to  be  presented  later  in  this  report 
indicate  that  any  such  errors  were  not  important.  The  "transition-end" 


criterion  used  for  the  various  records 


(iF)  <  2  (iF)upstream  and  (If)  <  l-5  (!f)upstr< 
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Sample  displays  of  velocities,  Reynolds  stresses  and  the  turbulence 
detector  function  are  presented  in  Figures  14  through  24.  These  particular 
samples  were  selected  to  display  the  performance  of  the  turbulence  detector 
for  a  wide  range  of  intermittency  levels  and  freestream  turbulence  intensi¬ 
ties.  In  order  to  present  the  time  resolution  of  events,  time  blocks  of  only 
200  ms  are  displayed.  These,  then,  represent  only  4%  of  the  5  seconds  of  con¬ 
tinuous  sample  obtained  for  each  total  record.  The  location,  acceleration 
level,  freestream  turbulence  intensity  and  intermittency  level  (over  the  en¬ 
tire  5  second  record)  are  given  for  each  figure. 


Sample  ver  ical-x-wire  and  single-wire  records  taken  at  Y/6  =  0.275  far  3 
low  freestream  turbulence  intensity,  moderate  int e rm  it t enc y  flow  are  presented 
in  Figure  14.  Note  that  the  x-wire  (14a)  and  the  single-wire  (14b-)  data  were 
not  obtained  simultaneously  so  the  u-t  ime  records  should  not  correlate. 

Figure  14a  gives  the  u,v  and  Reynolds  stress  records  and  the  detector  function 
based  on  the  uv  product.  Starts  of  turbulent  periods  correspond  to  the  + 
spikes  of  the  detector,  stops  to  the  -  spikes.  Figure  14b  gives  the  single¬ 
wire  u-time  record  and  the  detector  function  based  on  u.  For  both  14a  and  14b 
the  starts  and  stops  of  turbulence  as  selected  by  the  conditional  sampling 
routine  coincide  with  distinctly  active  portions  of  the  records.  The 
hor izontal-x-wire  data  taken  at  the  same  location  and  conditions  is  for  Figure 
14  are  presented  in  Figure  15  (Fig.  15b  repeats  Fig.  14b  for  convenience).  As 
can  be  seen  from  an  examination  of  Figure  15a,  the  uw  discriminator  al  so 
accurately  identifies  the  starts  and  stops  of  turbulent  periods  .if  the  flow. 

Figures  16  and  17  are  similar  to  Figures  14  and  15  in  that  they  present 
uv ,  uw,  and  u  records  for  a  single  location  in  a  test  boundary  layer.  In  this 
case,  however,  the  samples  show  the  performance  of  the  detectors  for  a  high 
intermittency  flow  (y  =  0.68).  Again,  all  the  detectors  ( uv ,  uw ,  and  u ' 
accurately  identify  periods  of  turbulence. 

Figures  18  through  23  are  an  abbreviated  version  of  the  turbulence- 

detector  displays  in  that  only  the  uv ,  uw,  and  u  records  are  presented. 

Figures  18  through  23  present  sample  records  for  a  wide  range  of  flow 

conditions;  Figure  18  —  low  T  and  low  y,  Figure  19  —  moderate  T  and  low  y, 

E  h 

Figure  20  —  high  T  and  moderate  y,  Figure  21  —  low  T  and  high  y,  Figure  22 

E  E 

—  moderate  T^.  and  high  y  and  Figure  23  —  high  T  and  high  y.  Finally  Figure 
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24  presents  sample  single-wire  records  for  a  single  boundary  layer  at  three 
distinctly  different  distances  (Y/6  =  0.056,  0.193  and  0.926)  from  the  wall. 
Profiles  of  the  turbulence  intensity  and  intermittency  for  this  particular 
boundary  layer  are  given  in  Figure  25.  An  examination  of  Figure  24  and  25 
reveals  that  the  u  detector  function  was  effectively  able  to  identify 
turbulent  periods  dispite  the  fact  that  the  average  turbulence  level  varied 
enormously  across  the  boundary  layer.  It  should  be  pointed  out,  however,  that 
the  Reynolds  stress  detector  functions  performed  unsat isfactoral ly  in  the 
outer  part  of  the  boundary  layers.  For  Y/6  <  0.5  the  uv ,  uw,  and  u  criteria 
would  consistently  agree  within  about  ±  2%  on  the  level  of  intermittency. 
Beyond  this  point,  however,  for  most  of  the  time  the  levels  of  uv  and  uw  in 
the  boundary  layers  exceeded  the  upstream,  near-wall,  pre-bursting  "reference" 
value  and  the  "end-of-turbulence"  criterion  failed.  A  detector  using  a 
threshold  level  based  on  local  values  might  offer  some  promise  for  future 
invest igat ion . 

In  summary,  for  all  three  detector  functions  (for  Y/6  <  0.5)  the  condi¬ 
tional  sampling  program  was  demonstrated  to  be  a  highly  versatile  and  accurate 
procedure  for  detecting  turbulent  periods  in  transitional  boundary  layer 
f  low . 


A  A  A 


l  >V mVuVaVmV  a  a  „V 


ViV.  .Vaa'A*’’  -  ‘n.’’  \  V-*  V 


3 


-V. 7>lVrCTLV\>7-V  -v 


V.  EXPERIMENTAL  RESULTS 


1.  STREAMWISE  DISTRIBUTION  OF  INTERMITTENCY 


The  "universal"  distribution  of  Dhawan  and  Narasimha  (Ref.  25)  is  widely 


accepted  as  an  accurate  description  of  the  development  of  near-wall  intermit- 


tency  in  zero  pressure  gradient  transitional  boundary  layers. 


Y  =  1  -  exp  [-  0.412  C  J 


5  -  (X  -  Xt)  A 


=  X  (y  =  0.75)  -  X  (y  =  0.25) 


For  experimentally  measured  intermittencv  distributions  the  beginning  of 


transition,  X^,  is  most  accurately  determined,  as  per  Narasimha  (Ref.  26),  by 


plot  t ing 


F(y)  =  [-  ln( 1  -  Y)] 


against  streamwise  location  and  finding  the  X  (y  =  0)  intercept.  As  dis¬ 


cussed  in  Reference  26,  in  this  format  zero  pressure  gradient  intermittencv 


distribution  data  fall  on  a  single  straight  line. 


The  distribution  of  the  near-wall  intermittencv  measured  for  the  four 


test  cases  of  the  present  program  are  presented  in  Figure  26.  These  near-wall 


intermittencv  values  represent  a  local  average  of  the  values  inferred  from  the 


uv,  uw,  and  u  detectors  as  described  in  CONDITIONAL  SAMPLING  TECHNIQUES, 


Because  these  independent  measurements  consistently  agreed  with  each  other 
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within  about  ±  2%  for  Y/ 6  <  0.5,  their  average  is  thought  to  be  a  highly  reli¬ 
able  measure  of  the  near-wall  value.  Figure  26  reveals  that  for  cases  K.20C1, 
K.20G2  and  K.75G2  there  was  a  distinct  "kink"  in  the  intermit tency  distribu¬ 
tions.  This  same  two-stage  behavior  for  transition  in  favorable  pressure 
gradients  was  recently  reported  by  Narasimha  (Ref.  27)  where  the  break  in  the 
intermittency  distributions  was  referred  to  as  a  "subtransition".  Note  that 
there  was  no  evidence  of  a  "subtransition"  for  the  K.75G3  case,  probably  a 
result  of  the  extremely  high  disturbance  level  of  that  test. 

The  intermittency  data  of  the  present  study  are  plotted  in  the 
"universal"  coordinates  in  Figure  27.  The  data  are  plotted  in  two  ways;  (1) 
the  origin  (X^)  was  taken  as  the  x  intercept  of  the  data  (open  symbols)  and 
(2)  was  determined  by  extrapolating  the  intermittency  data  for  X  > 
"subtransition"  to  F(  y)  =  0  (solid  symbols).  Neither  method  resulted  in 
agreement  with  the  "universal"  zero  pressure  gradient  distribution  over  the 
entire  length  of  transition.  By  using  the  X^  determined  from  extrapolation 
(solid  symbols),  however,  all  the  data  beyond  the  subtransition  fell  on  the 
universal  curve  (including  the  K.75G3  data  which  had  no  subtransition  point). 
The  extrapolation  procedure  failed  for  the  K.20G2  case  because  the  subtransi¬ 
tion  occurred  beyond  y  =  0.25,  the  lower  limit  in  the  definition  of  X.  The 
results  presented  in  Figures  26  and  27  as  well  as  data  of  Reference  27  indi¬ 
cate  that  it  is  possible  for  acceleration  to  substantially  inhibit  the  devel¬ 
opment  of  the  intermittency  until  some  critical  stage  (the  "subtransition")  is 
reached.  Following  the  "subtransition"  point  the  intermittency  distribution 
appears  to  follow  the  same  universal  scaling  low  observed  for  zero  pressure 
gradient  transitions. 


2.  BOUNDARY  LAYER  TURBULENCE  PROFILES 


Boundary  layer  mean  and  fluctuating  velocity  profile  data  were  obtained 
at  approximately  seven  streamwise  locations  for  each  of  the  four  experimental 
test  cases.  Single-horizontal  hot  wire  (u  component  only)  data  were  obtained 
at  all  upstream  stations  where  the  test  boundary  layer  had  not  grown  to  suffi¬ 
cient  thickness  to  permit  the  use  of  x-wire  probes.  For  the  downstream  4 
stations  of  each  test  both  single-wire  and  vertical  and  horizontal  x-wire  (u, 
v,  and  w  components)  data  were  recorded.  As  an  example,  the  profile  data 
obtained  for  one  of  the  test  cases  (K  =  0.75  x  10-6,  GRID  #2,  K.75G2)  are 
presented  in  Figures  28A  through  28K  and  will  be  discussed  here.  The  profile 
data  for  all  four  test  cases  are  provided  in  both  plotted  and  tabulated  form 
in  Appendix  A. 

The  figures  presented  in  this  section  and  in  Appendix  A  were  all  gener¬ 
ated  using  a  computer  based  plotting  system  (GOLDEN  GRAPHER).  For  consis¬ 
tency,  all  of  these  figure  formats  (sizes,  scales,  labels,  symbols,  etc.)  were 
kept  the  same  for  the  entire  data  set.  For  each  of  the  downstream  locations 
10  small  plots,  five  per  figure,  of  various  fluctuating  quantities  are  pre¬ 
sented.  Since  only  a  few  quantities  were  measured  for  the  upstream  profiles, 
the  locations  on  the  figures  assigned  for  the  remaining  plots  appear  as  blank 
areas  on  the  page. 

The  profile  data  (u  component  only)  obtained  at  the  most  upstream  sta¬ 
tion,  X  =  8.8  inches,  are  given  in  Figure  28A.  At  this  station  there  were  no 
turbulent  bursts  detected  and  only  the  mean  velocity  and  long  term 
("COMPOSITE")  average  U/U  are  presented.  Even  though  there  were  no  periods 


when  the  turbulence  criterion  was  satisfied  the  turbulence  intensity  (u'/U^) 
reached  nearly  7%  at  about  Y/S  =  0.3.  This  location  for  the  maximum  amplitude 
of  fluctuations  is  in  excellent  agreement  with  the  maximum-amplitude  locations 
found  from  both  theoretical  predictions  and  experimental  data  for  low- 
amplitude,  To  1 lmien-Schl icht ing  disturbances  in  zero  pressure  gradient 
boundary  layers  (Ref.  28).  This  in  spite  of  the  fact  that  this  test  profile 
was  highly  subcritical.  For  the  profile  of  Figure  28A  Re  =  230  while  the 
critical  Reynolds  number  for  linear  disturbance  amplification  for  this 
acceleration  level  (Ref.  28)  was  Reg  =  1580. 

In  addition  to  the  composite  turbulence  intensity  distribution,  Figure 
28A  also  presents  the  long-term  ("composite")  mean  velocity  distribution.  For 
reasons  related  to  the  clarity  of  presentation  of  the  downstream  profiles  only 
straight-line  segments  connecting  measured  data  (no  data  symbols)  for  U/U^.  > 
0.6  are  plotted.  These  hot  wire  mean  velocity  profile  data  were  in  excellent 
agreement  with  similar  profile  data  measured  with  a  miniature  Pitot  probe  as 
part  of  Reference  16. 

The  results  for  the  next  streamwise  station,  X  =  12.8  inches,  are  given 
in  Figure  28B  where,  again,  there  were  no  turbulent  periods  detected.  The 
maximum  boundary  layer  turbulence  intensity  had  grown  to  about  8%  and  was,  as 
with  Figure  28A,  located  at  about  Y/5  =  0.3.  Note  that  the  local  freestream 
turbulence  intensity  above  the  station  is  given  for  each  figure  and  that  T 

E 

fell  from  2.0  to  1.9%  between  the  first  two  stations. 

For  the  next  station  (X  =  16.8  inches.  Figure  28C)  the  conditional  sampl¬ 
ing  procedure  determined  that  the  criteria  for  turbulence  were  met  for  a 
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small  fraction  (about  2%)  of  the  time.  Three  additional  plots  are  presented 
for  Figure  28C  (1)  the  intermittencv  (y)  vs  Y/6,  (2)  the  int er-t urbu 1 ent -zone 
intensity  vs  Y/6  and  (3)  the  turbulent-zone  intensity  versus  Y/6.  The  turbu¬ 
lent  and  inter-turbulent-zone  plots  present  conditional  (zonal)  averages  of 
the  quantities  during  and  between  turbulent  periods.  The  term  inter-turbulent 
was  chosen  instead  of  laminar  because  of  the  highly  disturbed  character  of  the 
flow  between  turbulent  bursts.  For  the  test  cases  with  high  freestream  turbu¬ 
lence  there  were  large-amplitude,  relatively  low-frequency  fluctuations 
between  the  relatively  high-frequency  turbulent  fluctuations.  For  example  see 
the  time  displays  of  Figure  20. 

For  Figure  28C,  the  turbulence  intensity  during  the  turbulent  periods  is 
greater  in  magnitude  and  peaks  much  closer  to  the  wall  than  the  composite 
profile.  The  composite  and  inter-turbulent  profiles,  meanwhile  show  an 
increase  (from  8  to  9%)  in  peak  intensity  from  the  previous  station.  The 
intermittencv  profile  indicated  that  this  quantity  was  practically  constant 
across  the  boundary  layer  at  about  2%.  Two  additional  data  fairing  lines,  the 
turbulent  and  inter-turbulent  zone  average  velocity  ratios,  can  also  be  seen 
in  the  Mean  Velocity  plot.  All  three  mean  velocity  ratios  are  plotted  against 
Y/6  where  6  is  the  long-term  mean  boundary  layer  thickness.  As  will  be  demon¬ 
strated  later  in  this  report,  the  turbulent  periods  correspond  to  regions  of 
greatly  increased  boundary  layer  thickness.  Since  the  boundary  layer 
thickness  during  the  turbulent  period  is  much  greater  than  the  long-term 
boundary  layer  thickness  the  turbulent-zone  mean  velocity  decreases  relative 
to  the  composite  average.  Since  the  intermit tency  was  so  low  for  this  case 
the  composite  and  inter-turbulent  average  were  practically  identical. 
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Figures  28D  and  E  present  the  first  full  set  of  profile  data  including 
results  from  both  single-wire  and  x-wire  probes.  The  Intermittencv  distribu¬ 
tion  plot  of  Figure  28D  indicate  that  all  three  conditional  detectors  indi¬ 
cated  that  y  *  0.08  across  most  of  the  boundary  layer.  Sample  velocity  data 
records  for  this  profile  at  Y/S  =  0.286  can  be  seen  in  Figure  19. 

The  mean  velocity  and  streamwise  component  turbulence  intensity  results 
for  Figure  28D  were  very  similar  to  those  of  the  previous  station.  The  normal 
and  transverse  profiles  indicate  that  at  this  early  stage  of  transition  the 
turbulence  is  highly  anisotropic  with  a  particularly  small  v'  component.  Five 
additional  plots  must  be  introduced  for  Figure  28E.  The  Shear  Stress  and 
Turbulence  Kinetic  Energy  plots  display  the  distribution  of  the  conditional 
and  composite  averaged  values  of  those  quantities  across  the  boundary  layer. 
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Figures  28F/G,  28H/I  and  28J/K  present  a  similar  set  of  10  plots  for  X  = 
28.8,  36.8,  and  48.8  inches.  Examination  of  these  figures  reveals  that  with 
increasing  X  the  boundary  layer  turbulence  structure  progres s ive 1 v  approached 
that  expected  for  equilibrium  turbulent  flow.  The  various  intermittency 
detector  were  generally  in  very  good  agreement  for  Y/5  <  0.5.  As  previously 
discussed  the  uv  and  uw  detector  functions  performed  poorly  in  the  outer  por¬ 
tion  of  the  boundary  layers. 

Similar  sets  of  turbulence  profile  data  were  obtained  for  the  remaining 
three  test  combinations  of  acceleration  level  and  turbulence  intensity  (see 
APPENDIX  A).  In  an  attempt  to  reach  general  conclusions  about  the  development 
of  the  turbulence  through  transition  the  four  structural  coefficients  measured 
for  all  of  the  test  cases  were  combined  in  Figure  29.  In  this  figure  the 
structural  coefficients  are  plotted  as  a  function  of  the  local  int e rmi t t enc v 
at  that  particular  station.  The  individual  coefficients  are  normalized  bv  the 
value  measured  for  a  fully  turbulent  equilibrium  boundary  layer  at  the  same 
freestream  turbulence  intensity  in  the  study  of  Ref.  29.  Each  figure  includes 
the  normalized  coefficients  based  on  both  the  composite  and  turbulent-zone 
conditional  averages.  Although  considerable  scatter  is  evident  when  the  data 
are  plotted  in  this  format  it  is  possible  to  draw  three  conclusions  from  the 
trends.  (1)  For  all  four  plots  the  present  data  indicate  that  the  structural 
coefficents  approached  the  equilibrium  turbulent  boundary  laver  values 
(Ref.  29)  as  the  intermittency  reached  100  percent.  (2)  The  turbulent-zone 
conditionally  averaged  structural  coefficients  showed  much  less  variation  with 
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interraittency  than  did  the  composite  values.  In  fact  the  turbulent-zone  con¬ 
ditionally  averaged  direct  stress  coefficients  (a,  b,  and  c)  were  nearlv  equal 
to  the  fully  turbulent  values  for  all  intermi t tenc ies .  (3)  The  trends  of 

coefficients  a  and  b  reveal  the  extreme  anisotropy  of  the  turbulence, 
especially  in  the  early  stages  of  transition.  Commonly  used  equilibrium 
values  for  a  and  b  are  0.5  and  0.2  respectively  giving  a  ratio  a/b  =  2.5.  The 
present  data  indicate  that  at  the  early  stages  of  turbulent  bursting  this 
ratio  was  approximately  13. 


3.  CONDITIONALLY  AVERAGED  MEAN  VELOCITY  PROFILES 

Momentum  and  displacement  integral  thicknesses  were  computed  for  each  of 
the  conditionally  averaged  mean  velocity  profiles.  For  each  intermittent 
profile  integral  thicknesses  were  determined  for  the  turbulent  zone,  the 
inter-turbulent  zone  and  the  composite  average.  It  was  observed  that  the 
inter-turbulent  profiles  were  laminar-like  with  shape  factors  (H  =  6  / Z)  gen¬ 
erally  equal  to  2.0  or  greater.  Because  of  their  laminar-like  shapes  skin 
friction  coefficients  were  computed  for  these  inter-turbulent  zone  profiles 
using  the  molecular  viscosity  and  the  near-wall  normal  velocity  derivative. 

It  was  also  observed  that  the  turbulent  zone  profiles  were  turbulent-like  with 
shape  factors  equal  to  about  1.5.  When  plotted  in  the  universal  turbulent 
coordinates,  U  +  versus  Y  +,  these  turbulent  zone  profiles  could  be  well 
represented  by  the  turbulent  law-of-the-wal 1 .  Skin  friction  coefficients  for 
these  turbulent  zone  profiles  were  determined  by  a  least-squares  best  fit  to 
the  standard  log-wall  relationship 
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The  conditionally  averaged  mean  velocity  profile  results  for  the  four 

test  cases  are  presented  in  Figures  30  through  33.  For  each  figure  the  shape 

★ 

factor  6  /0,  momentum  thickness,  0,  Reynolds  number  based  on  0  and  skin  fric¬ 
tion  coefficient  distributions  are  plotted  as  a  function  of  streamwise  dis¬ 
tance  along  the  test  plate.  The  turbulent-zone,  inter-turbulent-zone  and 
composite  values  are  given  for  each  plot.  The  composite  skin  friction  coeffi¬ 
cients  were  computed  by  weighting  the  turbulent  and  inter-turbulent  contribu¬ 
tions  according  to  the  local  intermittency . 

Cf  =  yCf  +  (1-y)  Cf  (18) 

rcomp  rturb  rinter-turb 

The  solid  symbols  downstream  of  the  present  transitional  results  indicate 
data  obtained  for  these  same  test  conditions  using  miniature  Pitot  tubes 
(Ref.  16).  These  transitional  mean  profile  data  will  be  compared  with 
predictions  from  boundary  layer  computation  procedures  in  a  later  section  of 
this  report. 


4.  ENSEMBLE  AVERAGED  TURBULENT- ZONE  PROFILES 


Ensemble  averaging  techniques  were  emploved  to  examine  the  variation  of 
the  mean  streamwise  velocity  profiles  through  the  turbulent  periods  of  the 
flow.  An  ensemble  averaging  data  reduction  program  was  constructed  to 
uniquely  identify  each  individual  turbulent  period  and  determine  the  time 


interval  between  its  beginning  and  end. 


determined  by  the  burst  detector) 


The  "start"  of  a  turbulent  period  (as 
is  defined  as  t  =  0  for  these  ensemble 


computations  while  the  time  period  from  beginning  to  end  is  designated  as  At. 


For  each  detected  turbulent  period  the  velocity  versus  time  record  was  bro- 

At 

ken  into  intervals  of  yjj  with  an  average  velocity  computed  for  each  interval. 
For  each  burst  interval-averaged  velocities  were  computed  for  five  intervals 
prior  to  t  =  0,  ten  intervals  during  the  turbulent  period  and,  finally,  for 
five  intervals  following  t  =  At.  The  velocities  at  each  particular  fraction 
of  a  turbulent  period  averaged  over  all  the  turbulent  periods  constituted  the 
ensemble-averaged  local  velocity.  These  ensemble-averaged  velocities  were 
constructed  from  the  individual  velocity  versus  time  record  (5  seconds  long) 
representing  the  instantaneous  velocity  at  a  fixed  location  from  the  wall. 

The  ensemble  routine  computed  the  ensemble-averaged  mean  velocity  variation 
through  the  turbulent  periods  for  that  particular  distance  from  the  wall. 

A  single  probe  at  a  fixed  location  in  a  transitional  boundary  laver  de¬ 
tects  a  wide  range  of  turbulent  period  lengths  (Ax's).  These  differing  peri¬ 
ods  can  correspond  to  the  passage  of  turbulent  patches  of  different  size  or  to 
random  transverse  locations  of  the  patches  relative  to  the  probe  (the  edge  of 
a  large  patch  clipping  the  probe  appears  as  a  short  turbulent  period!.  During 
the  examination  of  these  ensemble  averaged  distributions  it  was  observed  that 
the  results  were  much  more  coherent  if  only  the  burst  lengths  within  about  the 
middle  60%  of  the  length  distribution  were  included  in  the  ensembles.  The 
very  short  and  very  long,  about  the  extreme  20%  of  each  end  of  the  distribu¬ 
tion,  turbulent  periods  were  excluded  from  the  computations.  For  most  of  the 
results  presented  here  about  50  to  100  turbulent  periods  were  included  in  each 
ensemble  distribution. 
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In  order  to  reveal  the  presence  or  lack  of  general  trends  in  these 
results,  turbulent-zone  mean  velocity  profile  variations  were  examined  for 
four  distinctly  different  test  conditions.  These  four  particular  cases  were 
selected  based  upon  two  criteria:  (1)  each  case  corresponded  to  a  different 
one  of  the  four  accelerat ion/turbulence  test  combinations  and  (2)  each  case 
corresponded  to  a  different  intermit tency  level.  The  ensemble-averaged  veloc¬ 
ity  distributions  for  these  four  cases  are  presented  in  Figures  34  through  37. 
The  following  Table  lists  the  test  conditions,  locations  and  intermittencies 
for  these  figures. 
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There  are  a  number  of  features  common  to  these  four  figures.  The  ordi¬ 
nate  for  Figures  34  through  37  is  the  time-interval  ensemble-averaged  velocitv 
(<  U  >)  normalized  by  the  freestream  mean  velocity.  For  each  profile,  veloc¬ 
ity  histories  are  given  for  20  normal  locations  ranging  from  about  Y/5  =  0.05 
to  about  Y/6  =  .83.  The  ensemble-averaged  velocities  are  plotted  as  a  func¬ 
tion  of  turbulent  period  fraction  (t/Ai).  Zero  corresponds  to  the  "start" 
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time  of  the  detector  function.  It  is  important  to  point  out  that  these 
ensemble-averaged  plots  are  based  on  locally  determined  timing  of  the  turbu¬ 
lent  periods.  The  beginning  and  end  of  turbulent  pe  iods  was  determined  at 
that  specific  distance  from  the  wall.  This  technique  does  not  consider  the 
shape  of  the  leading  and  trailing  faces  of  the  turbulent  patches  in  that  all 
the  t  =  0  and  u  =  At  events  were  plotted  as  if  they  were  simultaneous. 

The  most  distinct  feature  of  Figures  34  through  37  is  the  sharp  increase 
in  velocity  near  the  wall  at  t  =  0  followed  by  a  fall  to  about  the  upstream 
value.  The  velocity  histories  for  the  outer  portion  of  the  boundarv  lavers 
indicates  a  significant  deceleration  at  the  arrival  of  the  turbulent  patch 
followed  by  a  gradual  return  to  the  inter-turbulent  velocity.  Discussion  of 
the  individual  cases  follows. 

The  velocity-history  plots  of  Figure  34  are  for  a  case  with  a  very  high 

freestream  turbulence  intensity  (T  =  4.7%)  and  a  very  low  near-wall  intermit- 

E 

tency  level  (y  =  0.12).  The  ensemble-averaged  velocity  histories  of  Figure 
34,  then,  portray  very  small,  young  turbulent  patches  in  a  very  highly  dis¬ 
turbed  environment.  This  case  evidenced  the  lowest  degree  of  "organization" 
of  the  four  cases  studied.  Note  that  the  upstream  (t  <  0)  and  downstream  (t  > 
At)  velocities  do  not  agree  well  for  all  normal  locations.  The  apparent 
"disorganization"  of  the  flow  of  Figure  34  probably  indicates  that  the 
individual  turbulent  patches  were  still  developing  at  this  near-onset 
locat ion . 


Figure  35  presents  a  case  with  moderate  (T  =  1 . 7"0  freestream  turbulence 

and  high  intermittency  (y  =  0.78).  The  individual  velocity  histories  are  much 

smoother  and  more  correlated  for  Figure  35  than  for  Figure  34  indicating  that 

there  was  a  greater  burst -to-burst  consistency.  The  upstream  (t  <  0)  and 

downstream  (t  >  Ax)  velocities  are  not  constant  for  Figure  35,  a  consequence 

of  the  high  intermittency.  At  this  late  stage  of  transition  the  inter- 

turbulent  intervals  were  small  compared  to  Ax  (see  Figs.  17  and  21).  3v 

Ax 

including  samples  -  — —  upstream  and  downstream  of  the  edges  of  the  turbulent 
periods,  random  samples  from  neighboring  turbulent  patches  have  been  lumped 
into  the  inter-turbulent  ensembles. 

Figure  36  presents  results  quite  similar  to  those  of  Figure  35.  In  this 
case,  however,  the  intermittency  level  was  much  1  'w.-r  (  <  =  0.31''  and  the  up¬ 
stream  and  downstream  velocities  are  both  constant  and  in  near  agreement. 

Finally,  a  case  with  very  low  freestream  turbulence  (T  =  0.7-v;)  is  given 
in  Figure  37.  Again,  because  of  the  high  intermittency  (y  =  0.68),  "leakage" 
from  neighboring  turbulent  patches  has  contaminated  the  upstream  and  down¬ 
stream  ensembles.  Tlie  velocity  histories  of  Figure  37  were  cross-plotted  into 
"instantaneous"  ensemble-averaged  profiles  for  Figure  3S.  The  ensemble- aver¬ 
aged  ratios  are  plotted  against  Y/5  where  6  is  the  long-term  average  of  the 
boundary  layer  thickness.  The  relative  time  (in  burst-length  fraction,  I-'At) 
for  each  profile  is  given  near  the  boundary  layer  edge.  Note  the  enormous 
increase  in  boundary  layer  thickness  and  change  in  profile  shape  as  the  turbu¬ 
lent  patch  passes  by.  Selected  profiles  from  Figure  38  are  presented  in 


Figure  39  in  the  universal  turbulent  profile  coordinates.  The  upstream  and 
downstream  profiles  (t/Ax  =  -  0.3  and  1.3)  are  seen  xo  be  in  excellent 


agreement  with  each  other  and  to  exhibit  a  standard  laminar  shape.  At  t  =  0 
U  +  >  Y  +  for  a  few  of  the  near  wall  points  hinting  of  a  possible 
instantaneous  adverse  pressure  gradient.  By  t/At  =0.1  the  profile  begins  to 
appear  transitional  and  by  the  midpoint  (t/At  =  0.5)  it  agrees  very  well  with 
the  fully  turbulent  law-o f-the-wai 1 .  The  variation  of  selected  profile 
quantities  for  the  "instantaneous"  ensemble-averaged  profiles  are  plotted  in 
Figure  40.  Figure  40  presents  the  overall  boundary  layer  thickness,  6,  the 
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momentum  thickness,  0,  the  shape  factor,  6  /0  and,  finally  the  skin  friction 
coefficient.  As  indicated  in  the  skin  friction  plot  these  instantaneous  skin 
friction  values  were  computed  from  the  laminar  viscosity  for  the  inter- 
turbulent  profiles,  from  the  1 aw-o f-the-wal 1  for  the  turbulent  profiles  and  by 
a  simple  numerical  average  of  the  two  for  the  others.  The  results  presented 
in  Figure  40  reveal  the  extreme  discrepancies  between  the  inter-turbulent  flow 
and  that  in  the  turbulent  patches.  A  turbulent  patch  appears  to  be  typically 
3  times  thicker  than  the  surrounding  inter-turbulent  boundary  layer  with  a 
momentum  thickness  nearly  twice  as  large  as  the  interturbulent  value.  The 
shape  factor  is  seen  to  plunge  from  about  2.3  to  about  1.6  and  then  return  to 
the  upstream  value.  The  skin  friction  is  nearly  5  times  greater  at  midpoint 
in  the  turbulent  patch  than  computed  for  the  inter-turbulent  profile.  Note 
that  these  "instantaneous"  ensemble-averaged  profile  results  are  in  very  good 
agreement  with  the  average  conditional  results  (T-turbulent ,  IT-inter- 
turbulent,  C-composite)  marked  on  the  plots. 

5.  SPECTRAL  DISTRIBUTIONS 

A  Fast-Four ier-Transform  (FFT)  routine  contained  in  the  Cranfield  Data 
Systems  package  was  employed  to  compute  the  spectral  distribution  of 
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individual  ve loc i t y- t ime  data  records.  The  previously  described  conditional 
sampling  program  was  configured  such  that  spectral  distributions  could  be 
computed  for  the  total  (composite)  data  sample  or  for  just  the  turbulent  or 
inter-turbulent  portions  of  the  record.  The  individual  inter-turbulent  and 
turbulent  periods  were  arranged  end-to-end  to  form  constructed  continuous 
records.  Although  this  procedure  does  not  compromise  resolution  of  the  high- 
frequency  portions  (upper  limit  =  sampling  rate/2)  of  the  spectrum  it  does 
restrict  the  lower  frequency  limit.  Only  for  the  continuous  records 
(composite)  can  it  be  assured  that  the  lowest  frequency  resolved  is  equal  to 
the  reciprocal  of  the  total  sample  length.  Since  all  sampling  for  this  pro¬ 
gram  was  conducted  at  7813  Hz/channel  the  upper  limit  of  frequency  resolution 
for  all  cases  was  3906  Hz. 

In  order  to  demonstrate  the  performance  of  the  spectral  analysis  svstem 
used  here  the  power-spectral-density  (PSD)  distribution  computed  for  a  post- 
transitional  (fully  turbulent)  boundary  layer  is  given  in  Figure  41.  This 
figure  gives  the  PSD  for  the  streamwise  fluctuation  recorded  near  the  wall  at 
a  station  downstream  of  transition  (y  =  1.00)  for  one  of  the  test 
configurations.  For  this  figure  the  spectral  distribution  is  plotted  in  non- 
dimensional,  turbulent  boundary  layer  co-ordinates  as  suggested  by  Bradshaw 
(Ref.  30).  The  agreement  between  the  present  data  and  the  results  of 
Reference  30  are  seen  to  be  excellent. 

The  primary  objective  of  the  present  spectral  measurements  was  to  obtain 
information  about  the  nature  of  the  disturbance  distribution  in  transitional 
boundary  layers.  In  order  to  achieve  this,  spectral  distribution  measurements 
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were  obtained  for  boundary  layers  at  various  levels  of  intermittencv  as  well 
as  for  upstream  pre-trans it ional  profiles  in  which  no  turbulent  periods  were 
detected.  These  present  transitional  boundary  layer  spectral  distribution 
data  provide  an  opportunity  to  compare  the  predictions  of  linear  stability/ 
amplification  theory  with  experimental  measurements.  It  should  be  pointed 
out,  however,  that  the  experimental  test  cases  studied  for  this  program  all 
involved  relatively  high  levels  of  freestream  turbulence  with  accelerating 
flow.  In  contrast,  linear  small-disturbance  amplification  theory  and  experi¬ 
ments  (Refs.  31  and  32)  have  focused  primarily  in  single  frequency  distur¬ 
bances  to  boundary  layers  with  very  low  freestrean  turbulence.  These  present 
spectral  data/linear  theory  comparisons  are  particularly  interesting  because 
the  present  test  cases  fall  into  the  important  category  of  "bypass  mode"  tran¬ 
sition  (Ref.  33) . 

Figures  42  through  49  present  spectral  distributions  measured  at  succes¬ 
sively  increasing  strearawise  locations  through  the  transition  process  for  the 
test  configuration  K  =  0.2  x  10'6,  Grid  No.  1.  For  these  figures  the  spectral 
distributions  are  presented  in  dimensional  coordinates  since  the  appropriate 
length  scale  is  unclear.  Figure  42  presents  the  spectral  data  for  the  most 
upstream  profile  (X  =  12.8  in)  at  which  no  turbulent  periods  were  detected. 

At  this  station  the  momentum  thickness  Reynolds  number,  Re^,  was  360,  far 
below  the  critical  (Re^  =  680)  thickness  for  this  acceleration  level 
(Ref.  34).  All  frequencies  are  stable  at  this  Re^for  sufficiently  small  dis¬ 
turbances.  Spect-al  distributions  obtained  at  four  distances  from  the  wall 
are  given  in  Figure  42.  The  freestream  spectrum,  given  here  for  reference,  is 
in  excellent  agreement  with  the  distribution  measured  previously  by  Sato  for 
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this  turbulence  generating  grid  Reynolds  number  (Ref.  35).  The  three  spectra 
measured  in  the  profile  show  a  much  faster  roll-off  with  increasing  wavenumber 
than  expected  for  developed  turbulence.  There  was  a  ratio  of  nearly  6  orders 
of  magnitude  between  the  low  and  high-frequency  contributions  to  the  spectra 
at  both  Y/6  =  0.14  and  0.18.  These  results  indicate  that  deep  in  the  boundary 
layer  the  spectra  are  radically  different  from  the  broadband,  developed  free- 
stream  distribution  (the  driving  disturbance).  There  is  no  evidence  of  nar¬ 
rowband  disturbance  frequency  preference  at  this  profile.  Note  that  at  the 
two  locations  closest  to  the  wall  a  spike  corresponding  to  the  tunnel  blower 
blade-passing-frequency  can  be  discerned  (BPF).  This  spike  is  almost  certain¬ 
ly  an  acoustic  wave  and  not  a  velocity  fluctuation. 

Spectral  distributions  for  the  next  station  downstream  (X  =  16.8  in)  are 
given  in  Figure  43.  Again  the  freestream  spectrum  is  included  for  reference. 
The  near-wall  spectrum  has  changed  radically  from  the  previous  station  with 
four  bands  of  highly  amplified  disturbances.  The  band  labled  FI  corresponds 
to  the  most-unstable  frequency  predicted  by  linear  theory  for  this  accelera¬ 
tion  level  (Ref.  28).  The  2nd,  3rd,  and  4th  harmonics  of  FI  also  match  the 
observed  amplified  bands.  This  agreement  is  rather  remarkable  in  that  the 
boundary  layer  is  still  highly  sub-critical.  The  observed  amplified  frequen¬ 
cies  corresponded  to  the  predicted  critical  frequencies  for  this  large- 
amplitude  disturbance  level  even  though  the  linear  theory  would  predict  that 
Re^  was  subcritical  for  all  wavenumbers.  Note  that  the  BPF  acoustic  wave  is 
still  discernable. 

At  the  next  station  (X  =  20.8  in,  Fig.  44)  the  near-wall  disturbances 
have  become  both  more  powerful  and  broadband.  There  is  still  some  meager 


evidence  of  preferred  amplification  for  FI  and  2F 1 .  Note  that  at  this  station 
some  turbulent  periods  were  detected  with  y  =  0.002.  The  same  trend  continued 
for  Figure  45  with  a  steadily  increasing  contribution  from  the  higher  wave- 
numbers  to  the  near-wall  spectrum.  At  this  station  y  =  0.013. 

The  format  for  the  next  four  figures  (46  through  49)  differs  from  that  of 
Figures  42  through  45.  For  these  four  downstream  stations  the  spectra  from 
only  a  single  normal  location  (Y/S  =  0.10)  are  given  for  each  streamwise  sta¬ 
tion.  For  these  downstream  stations  the  intermittency  level  was  high  enough 
so  that  adequate  conditionally  sampled  turbulent  spectra  could  be  computed. 

For  these  four  profiles,  then,  spectra  are  given  for  the  composite  record  and 
for  the  inter-t urbulent  and  turbulent  portions  of  the  record.  The  axes  coor¬ 
dinates  are  subscripted  with  i  corresponding  to  either  the  composite,  inter- 
turbulent,  or  turbulent  value.  The  conditionally  sampled  spectra  show  the 
distributions  of  the  fluctuations  within  that  particular  zone.  The  turbulent 
spectra  for  all  four  streamwise  locations  appear  to  indicate  fully  developed 
turbulence.  For  example,  the  turbulent  spectral  distribution  is  seen  to  agree 
well  with  a  -5/3  slope  for  Figure  48.  The  inter-turbulent  spectra  change 
dramatically  with  increasing  streamwise  location.  At  X  =  28.8  in.  (Fig.  46) 
there  is  a  ratio  of  about  5  orders  of  magnitude  between  the  low  and  high  fre¬ 
quency  contributions.  By  X  =  52.8  in.  (Fig.  49)  the  relative  power  of  the 
high  frequency  content  has  greatly  increased  so  that  this  ratio  has  dropped  to 
about  3  orders  of  magnitude.  None  of  the  composite,  inter-turbulent  or  turbu¬ 
lent  spectra  show  any  evidence  of  narrow-band  preferred  amplification. 


The  final  two  spectral  distributions  presented  here  were  obtain*’ 1  for  the 
test  case  K  =  0.75  x  10-t3,  GRID  No.  2.  Figure  50  shows  distributions  obtained 
at  an  upstream  station  (X  =  8.8  in.)  where  there  were  no  detected  turbulent 
periods.  As  with  Figure  42  the  near-wall  spectrum  contained  verv  little  con¬ 
tribution  from  the  high-frequency  fluctuations.  Again  the  BPF  acoustic  wave 
stand  out  clearly.  At  X  =  12.8  in.  (Fig.  51)  the  near-wall  spectrum  contains 
much  more  power  in  the  high  frequency  range  with  some  evidence  of  selective 
amplification.  The  most-unstable  frequency  for  the  acceleration  associated 
with  the  data  of  Figure  51  (FI)  appears  to  roughly  coincide  with  one  of  the 
preferred  bandwidths. 


VI.  ASSESSMENT  OF  ANALYTICAL  PREDICTIONS  OK  TRANSITION 


Numerical  experiments  were  undertaken  to  assess  the  ability  of  current 

methods  for  predicting  heat  transfer  during  transition  in  accelerating  flows 

with  various  levels  of  freestream  turbulence.  The  analytical  predictions  were 

made  using  the  UTRC  ABLE  boundary  layer  code  developed  by  Edwards  et  al . 

(Ref.  34).  These  predictions  are  compared  with  the  experimental  data  base 

obtained  by  Blair  (Ref.  16)  for  two  different  rates  of  acceleration  and 

several  different  freestream  turbulence  levels.  Analytical  predictions  of 

transition  use  essentially  three  different  methods.  The  first  and  earliest 

(Refs.  13,  25,  35,  36,  37,  and  38)  use  empirical  correlations  for  the  start 

and  end  of  transition  usually  using  Reynolds  number  Re  based  on  boundary 
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layer  momentum  thickness  and  then  sora»  ir.Lermittency  distribution  to  produce  a 
smooth  transition  from  laminar  to  turbulent  flow.  Of  these  methods,  only  the 
first  two  predict  both  the  beginning  and  end  of  transition.  The  others 
predict  only  the  start  of  transition  (Ref.  36),  the  end  of  transition 
(Ref.  37),  or  the  intermittency  (Refs.  25  and  38). 

These  empirical  correlations  generally  suffer  from  difficulty  in  applica¬ 
tion  and  interpretation  in  flows  which  have  varying  freestream  conditions 
(i.e.,  flow  history).  Indeed  the  transition  from  laminar  to  turbulent  flow 
will  depend  on  flow  history  and  thus  these  methods  will  always  have  this 
problem.  In  view  of  this,  a  second  method  has  been  developed  which  depends  on 
the  solution  of  a  differential  equation  (usually  the  turbulence  kinetic  energy 
equation).  The  most  successful  of  these  was  developed  by  McDonald  and  Fish 
(Ref.  10)  and  improvements  were  later  made  by  McDonald  and  Kreskovsky 
(Ref.  39).  This  method  attacks  the  problem  of  transition  on  a  more 
fundamental  level  and  hence  presumably  is  applicable  to  a  wider  range  of 


problems  if  the  turbulence  has  been  properly  modeled.  However,  these  methods 
require  a  deeper  understanding  of  the  structure  of  turbulence  during  and  after 
transition  and  considerably  more  data  to  evaluate  empirical  constants.  These 
methods,  however,  have  an  advantage  over  methods  based  on  empirical  correla¬ 
tions  because  there  is  less  ambiguity  in  application  to  flows  with  a  varying 
flow  history  (i.e.,  varying  freestream  conditions  —  pressure  gradient  and 
freestream  turbulence  and  varying  wall  conditions  —  wall  temperature  or  heat 
flux) . 

With  the  development  of  the  two  equation  models  of  turbulence  (Refs.  40 
and  41)  came  the  possibility  of  predicting  transition  directly  (Refs.  42-47), 
since  these  models  solve  partial  differential  equations  for  the  turbulence 
properties.  Of  these  methods,  References  45  and  48  appear  the  best  documented 
and  appear  to  predict  many  of  the  qualitative  and  quantitative  features  of 
transition  in  both  constant  and  accelerating  freestream  flows.  An  attractive 
feature  of  these  methods  is  that  since  the  history  of  turbulence  is  modeled, 
transition  is  allowed  to  take  place  in  a  more  natural  manner.  The  accuracy  of 
such  predictions  for  general  flow  situations  depends  very  much  on  the 
universality  of  the  empirical  coefficients  used  in  the  various  models.  Thus, 
although  it  is  quite  possible  to  make  ad  hoc  adjustments  to  these  empirical 
constants  so  that  transition  can  be  predicted  for  a  limited  number  of  cases, 
it  is  not  clear  that  such  ad  hoc  changes  represent  a  phenomenological  modeling 
of  the  turbulent  processes,  particularly  when  nominal  constants  are  permitted 
to  change.  The  only  true  procedure  for  verifying  these  constants  is  to  obtain 
detailed  experimental  data  of  the  turbulent  parameters  so  that  all  empirical 
constants  can  be  phenomenologically  modeled.  At  a  minimum,  this  would  require 
measurements  of  dissipation,  a  difficult  but  not  impossible  task  so  that  both 
dependent  variables  in  the  <  -  e  model  are  known  experimentally.  Thus,  the 


third  method  for  predicting  transition  requires  an  even  deeper  understanding 

of  the  structure  of  turbulence  during  and  after  transition  and  even  more 

experimental  data  to  evaluate  empirical  constants. 

The  following  study  concentrates  on  only  the  first  two  methods  and  will 

be  limited  to  the  empirical  method  of  AbuGannom  and  Shaw  (Ref.  35)  (AGS)  and 

the  modified  one  equation  method  of  McDonald  and  Kreskovsky  (Ref.  39)  (MFK) . 

The  (AGS)  method  is  primarily  a  method  for  predicting  the  interm  it t ency 

distribution.  Since  this  report  presents  measurements  of  intermittency  y 

during  transition,  a  more  detailed  evaluation  of  this  method  can  be  made  than 

was  previously  possible.  The  (MKF)  method  is  primarily  a  method  for 

predicting  a  mixing  length  scale  Z  during  transition.  The  primary  empirical 

00 

factors  are  the  structual  coefficients  a,  b,  d  (especially  d)  which  vary 
during  transition.  Again  this  report  presents  measurements  of  these 
structural  coefficients  and  a  more  detailed  evaluation  of  this  method  is 
possible.  A  proper  evaluation  of  a  two  equation  <  -  e  model  of  turbulence  for 
predicting  transition  requires  at  a  minimum  the  measurement  of  dissipation  e. 
Unfortunately  this  data  is  not  available  at  the  present  time  and  so  this  model 
will  not  be  considered. 

1.  COMPARISON  OF  CALCULATIONS  WITH  EXPERIMENTAL  DATA 

The  cases  shown  in  the  following  figures,  are  all  presented  in  the  same 
manner.  In  addition  to  the  prediction  of  the  heat  transfer  during  transition 
from  laminar  to  turbulent  flow,  these  figures  also  contain  the  prediction  for 
completely  laminar  flow  (i.e.,  no  grid)  as  well  as  completely  turbulent  flow 
which  has  been  tripped  at  the  leading  edge.  In  all  cases  it  was  found  that 
the  thermal  history  of  the  boundary  layer  (i.e.,  wall  temperature  history)  was 


important.  Therefore  all  cases  were  calculated  using  the  measured  wall 
temperature  distribution.  For  the  completely  turbulent  boundary  layer  case, 
measurements  were  not  made  of  either  heat  transfer  or  wall  temperature,  there 
fore  the  wall  temperature  for  the  completely  laminar  flow  (no  grid)  case  was 
used.  This  latter  calculation  should  be  considered  as  a  reference  only.  The 
test  matrix  of  flow  conditions  to  be  examined  is  shown  in  Table  1.  They 
consist  of  two  different  pressure  gradients  K  =  0.20E-06  and  K  =  0.75E-06 


where 


K  - 

ir  dx 


and  four  different  turbulence  levels  (including  no  freestream  turbulence). 
The  nominal  upstream  turbulence  levels  are  given  on  Table  1  for  a  position 
12.0  in  upstream  of  the  leading  edge  where 


-  u'u' 


The  matrix  of  test  conditions  calculated  is  shown  on  Table  2  showing  the 

combinations  of  accelerations  and  freestream  turbulence  levels.  The  results 

are  presented  in  terms  of  the  variations  of  Stanton  number  (St  )  with  X 

e 

where 


ste  = 


pUeCp(Tw  -  V 


It  is  noted  that  wall  is  not  heated  for  the  first  2.19  in  and  therefore 
adiabatic  wall  is  assumed. 


2 .  ABOGANNOM-SHAW  PREDICTIONS 

Figure  52  shows  a  comparison  of  the  predictions  using  the  (AGS)  model 
with  experimental  data  for  both  the  weak  and  strong  acceleration  cases.  For 


Table  1:  Matrix  of  Measured  Flow  Conditions 


Weak  Acceleration  Strong  Acceleration 
K  =  0.20E0  -  6  K  =  0.75E  -  06 


No  Grid 


Grid  1 


Grid  2 


Grid  3 


0.000 


0.012 


0.026 


0.000 


0.027 


0.066 


Table  2:  Matrix  of  Calculated  Flow  Conditions 


K  =  0.20E  -  06  K  =  0.75E  -  06 


d  2 
Grid  3 
laminar 
turbulent 


Grid  2 
Grid  3 
laminar 
turbulent 


Grid  2 
Grid  3 
1 aminar 
turbulent 


Grid  2 
Grid  3 
laminar 
turbulent 


(AGS)  AbuGannom-Shaw  [3]  > 

Grid  1 

Grid  2 
laminar 

turbulent 

(MFK)  Mcdonald-Kreskovsky  [10] 

Grid  1 

Grid  2 

laminar 

turbulent 

(INP)  Experimental  Intermi t tency 

Grid  1 

Grid  2 

1 aminar 

turbulent 

(INA)  Experimental  Structural  Coefs. 

Grid  1 

Grid  2 

1  ami nar 

turbulent 

these  cases,  the  Cebec i-Smith  (Ref.  49)  eddy  viscosity  turbulence  model  was 
used.  Since  the  flow  has  varying  freestream  conditions,  the  following  proce¬ 
dure  was  used  to  predict  transition.  Let  Re^  be  the  local  Reynolds  number 
based  on  momentum  thickness  and  let  Re„  be  the  predicted  start  of  transition 
using  theAbuGannom-Shaw  transition  model  and  local  conditions.  Transition  was 
assumed  to  occur  when  Re^  =  Re^.  It  can  be  seen  from  these  figures,  that  the 
(AGS)  transition  model  produces  a  fairly  good  prediction  for  the  start  of 
transition,  but  poorly  predicts  the  end  of  transition.  Clearly  the 
intermittency  distribution  is  not  predicted  well  for  these  cases. 


3.  MCDONALD  -  FISH  -  KRESKOVSKY  PREDICTIONS 


Figure  53  shows  a  comparison  of  the  predictions  using  the  (MFK)  model 
with  experimental  data  for  both  the  weak  and  strong  acceleration  cases.  For 
these  cases,  the  experimental  freestream  turbulence  was  also  input.  As  can  be 
seen  from  these  comparisons,  the  (MFK)  model  does  not  predict  the  start  of 
transition  nor  does  the  model  fully  capture  the  fully  turbulent  flow  region. 

A  comparison  of  the  results  for  the  weak  acceleration  Figure  53a  with  the 
results  for  the  strong  acceleration  Figure  53b  appear  to  indicate  that  this 
model  overestimates  the  impact  of  acceleration.  It  should  be  noted  that  at 
the  time  that  the  model  was  developed,  the  structural  coefficients  a,  b,  d, 


given  by  equations  12,  13,  and  15,  and  the  dissipation  length  scale  L  given 

d 


L£j  =  (-u'u')  /  e 


were  not  known  for  transitional  flow,  and  hence  various  assumptions  were 


made  to  obtain  them  by  comparison  of  mean  flow  predictions  with  experimental 
data . 


4b 


4.  CEBECI-SMITH  PREDICTIONS  WITH  EXPERIMENTAL  INTERMITTENCY 


Since  the  intermittency  distribution  was  measured  for  these  cases,  the 
flow  field  and  heat  transfer  were  recalculated  using  the  measured  intermit¬ 
tency  distributions  together  with  the  Cebeci-Smith  (Ref.  49)  turbulence  model 
The  results  are  shown  on  Figures  54  through  57  and  the  predictions  for  heat 
transfer,  skin  friction,  and  momentum  thickness  appear  to  be  quite  good.  One 
may  then  conclude  from  these  results  that  the  use  of  the  intermittency 
distribution  is  a  valid  method  for  providing  a  smooth  transition  from  laminar 
to  turbulent  flow  and  that  if  one  had  good  correlations  for  the  beginning  and 
end  of  transition  quite  good  predictions  could  be  made  even  though  no 
provision  is  made  in  the  Cebeci-Smith  turbulence  model  for  the  effects  for 
freestream  turbulence. 

5.  MCDONALD-KRESKOVSKY  PREDICTIONS  WITH  STRDCTURAL  COEFFICIENTS 

Empirical  correlations  for  the  structural  coefficients  a,  b,  d  as  a  fjnc 
t ion  of  intermittency  y  are  given  by; 

a  =  0.600  [  1.30  +  (  1.00  -  1 . 30  ) y]  (23) 

b  =  0.250  [0.25  +  (1.00  -  0 . 25 ) y ]  (24) 

d  =  0 . 1 2 1 5y 1 /2  (25) 
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These  correlations  were  programmed  in  the  ARLF  code  and  the  cases  rerun 
with  the  experimental  intermittency  distributions.  The  predictions  for  heat 
transfer  are  shown  on  Figure  57.  These  predictions  show  some  improvement  over 
the  predictions  using  the  (MFK)  model  alone  and  give  encouragement  to  the 


basic  concept.  It  should  be  uoled  that  measurements  of  dissipation  length  L 


have  not  been  made  so  that  it  is  not  possible  to  completely  test  the 
procedure.  Since  an  improved  model  is  desireable,  it  is  stronglv  recommendei 
that  measurements  of  dissipation  e  be  made. 


Li 


VII.  CONCLUSIONS 


V, 

V, 

V 
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A  combined  experimental  and  analytical  program  has  been  conducted  to 
examine  transitional,  accelerating  boundary  layer  flows  with  high  levels  of 
freestream  turbulence.  This  program  was  designed  to  complement  a  parallel 
study  conducted  previously  under  AFOSR  funding.  The  earlier  program  focused 
on  measurement  of  transitional  heat  transfer  distributions  for  four  combina¬ 
tions  of  streamwise  acceleration  and  freestream  turbulence.  The  present  pro¬ 
gram  was  designed  to  document  the  boundary  layer  turbulence  structure  for  the 
same  four  test  conditions.  In  addition  spectral  distribution  data  were 
obtained  both  within  the  test  boundary  layers  and  in  the  freestream  flow.  The 
measurements  from  the  earlier  program  combined  with  the  present  results  pro¬ 
vide  a  comprehensive  examination  of  accelerating  transitional  boundary  layer 
flow.  Analytical  results  include  comparisions  between  these  various  measure¬ 
ments  and  a  variety  of  numerical  transitional  boundary  layer  computation  pro¬ 
cedures.  Specific  conclusions  reached  from  the  experimental  and  analytical 
programs  were  as  follows. 

1.  EXPERIMENTAL  PROGRAM 

1)  Documentation  of  the  flow  conditions  for  each  of  the  four  test  cases 
consisted  of  measurements  of  the  streamwise  freestream  velocity  (static 
pressure)  distributions,  all  three  components  of  the  turbulence  beyond  the 
edge  of  the  test  boundary  layers,  distributions  of  the  streamwise  integral 
length  scale,  the  spectral  distribution  of  the  freestream  turbulence  and, 
finally,  the  location  and  extent  of  the  boundary  layer  transition  zones. 


Measurements  of  all  these  quantities  were  obtained  for  both  the  previous  and 
present  experiments  (using  different  techniques  and  in  differing  detail!  . 
Agreement  between  the  present  and  earlier  data  was  excellent  assuring  that  the 
present  test  conditions  duplicated  those  of  the  earlier  AFOSR  contract. 

2)  The  hot-wire  conditional  sampling  routine  developed  as  part  of  this 
study  has  been  shown  to  be  consistently  capable  of  discriminating  between 
turbulent  and  non-turbulent  periods  in  transitional  boundary  layer  flow.  The 
discrimination  routine  proved  to  be  reliable  for  a  wide  range  of  freestream 
turbulence  intensities  and  interm  it tency  levels.  The  conditional  sampling 
routine  made  possible  computation  of  zonal  (turbulent  or  "inter-turbulent") 
mean  and  fluctuating  boundary  layer  flow  quantities  and  ensemble-averaged 
velocity  distributions  within  and  between  the  turbulent  patches. 

3)  The  streamwise  distributions  of  near-wall  intermittencv  measured  for 
the  various  test  cases  indicate  that  there  are  two  distinctly  separate  stages 
of  the  transition  process  for  accelerating  flow.  These  data  indicate  that  it 
is  possible  for  acceleration  to  substantially  inhibit  the  development  of 
interra  it tenc y  until  a  critical  point  (the  "subtransition")  is  reached. 
Following  the  "subtransition"  point  the  intermittency  distribution  appears  to 
follow  the  same  universal  scaling  law  observed  for  zero  pressure  gradient 
transitions.  It  appears  that  the  feasibility  of  accurate  computation  of 
transitional  boundary  layer  development  could  be  significantly  improved  if  a 
generalized  procedure  were  available  for  the  prediction  of  the  location  and 
extent  of  both  stages  of  this  two-step  process. 

4)  Profiles  of  the  conditionally  sampled  direct  and  Reynolds  stresses 
across  the  test  boundary  layers  were  obtained  at  several  streamwise  stations 


for  all  four  flow  conditions.  At  stations  upstream  of  intermittent  turbulent 
bursting  the  amplitude  distribution  of  the  streamwise  fluctuations  was  very 
similar  to  that  predicted  by  linear  small-disturbance  theory.  The  profiles 
obtained  at  locations  with  very  low  intermittencv  indicate  that  at  this  early 
stage  of  transition  the  turbulence  is  highly  anisotropic  with  a  particularly 
small  v'  component.  The  ratio  between  the  streamwise  and  normal  velocity 
fluctuations  was  over  five  times  larger  for  these  low  intermittency  flows  than 
for  equilibrium  turbulent  boundary  layers.  The  turbulence  components  steadily 
approached  the  equilibirium  turbulent  boundary  layer  arrangement  as  the  inter¬ 
mittency  approached  unity.  Suggested  functional  relationships  are  given  for 
the  distributions  of  the  turbulence  structural  coefficients  through  tran¬ 
sition. 

5)  Momentum  and  displacement  integral  thicknesses  and  skin  friction 
coefficients  were  computed  for  each  of  the  conditionally  sampled  mean  velocity 
profiles.  It  was  observed  that  the  turbulent  zone  profiles  could  be  well 
represented  by  the  f ul ly-turbulent  law-of-the-wal 1 .  It  was  also  observed  that 
the  inter-turbulent  zone  profiles  were  laminar-like  with  shape  factors  equal 
to  2.0  or  greater. 

6)  The  passage  of  "typical"  turbulent  patches  was  examined  by  computing 
relative  time  (in  burst-length  fractions),  ensemble-averaged  "instantaneous" 
mean  velocity  profiles.  In  this  format  the  inter-turbulent  zone  profiles  were 
shown  to  exhibit  the  standard  laminar  shape  while  the  mid-patch  turbulent  pro¬ 
files  agreed  very  well  with  the  turbulent  law-of-the-wal 1 .  These  results 
revealed  the  extreme  discrepancies  between  the  inter-turbulent  flow  and  that 


bl 


in  the  turbulent  patches.  A  turbulent  patch  appears  to  be  typically  three 
times  thicker  than  the  surrounding  inter-turbulent  boundary  layer  with  a 
momentum  thickness  nearly  twice  as  large  as  the  inter-turbulent  value.  The 
skin  friction  is  nearly  five  times  greater  at  midpoint  in  the  turbulent  patch 
than  computed  for  the  inter-turbulent  profile. 

7)  Spectral  analysis  of  the  fluctuations  in  a  boundary  layer  upstream  of 
the  onset  of  intermittency  have  revealed  a  band  of  highly  amplified  distur¬ 
bances  near  the  most-unstable  frequency  predicted  by  linear  theory.  This 
result  is  rather  remarkable  in  that  the  boundary  layer  at  that  location  was 
highly  subcritical.  The  result  reveals  the  preferred  amplification  of 
selected  strong  (non-linear)  disturbances  and  demonstrates  at  least  one  pos¬ 
sible  path  of  a  "bypass"  transition.  Conditionally  averaged  spectral  distri¬ 
butions  indicate  that  the  turbulent  zone  spectra  are  very  similar  to  those  of 
an  equilibrium  turbulent  boundary  layer  for  the  entire  length  of  transition. 


The  inter-turbulent  zone  spectra,  however,  exhibit  a  much  lower  fraction  of 
high  frequency  power  for  low  intermittenc ies  and  only  assume  an  equilibrium 
turbulent  boundary  layer  distribution  as  the  intermittency  approaches  unity. 


2.  ANALYTICAL  PROGRAM 


1)  An  assessment  of  the  AbuGannom-Shaw  (Ref.  35)  empirical  correlations 
for  the  prediction  of  transition  in  accelerating  flows  shows  that  this 
procedure  predicts  the  onset  of  transition  quite  well  but  does  not  predict  the 
length  of  transition. 


*  x VW  '1  V* 


2)  An  assessment  of  the  McDonald-Kreskovsky  (Ref.  39)  one  equation  model 


for  the  prediction  of  transition  in  accelerating  flows  shows  that  this 


procedure  does  not  predict  the  onset  of  transition  nor  does  it  fully  capture 


the  fully  turbulent  flow. 


3)  Use  of  the  experimental  intermittency  distribution  in  place  of 


empirical  correlations  produced  very  good  predictions  of  transition.  This 


result  indicates  that  the  modeling  of  transition  using  case-specific 


intermittency  distributions  is  a  valid  concept  and  that  better  empirical 


correlations  for  the  intermittency  could  produce  accurate  predictions  of 


trans it  ion . 


4)  Use  of  experimental  values  in  place  of  fully  turbulent  values  for  the 


structural  coefficients  in  the  McDonald-Kreskovsky  (Ref.  39)  one  equation 


model  produced  better  predictions  for  transition.  This  result  indicates  that 


better  correlations  for  the  structural  coefficients  and  the  dissipation  length 


could  perhaps  produce  better  predictions. 


v/.-.V.V.v.'.-.v; 


K'l 


K 


I 

I 

Ski*. 


VIII.  LIST  OF  SYMBOLS 


a  ,b  ,c  ,d 
cf 

CP 

K 


H 


L 

q 


d 
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Re 


Re, 


t 

T 

T 

v 

T 


u,v,w 


u'  ,v'  ,w' 


U 

+ 

u 

UT 

<u> 


u,2(f) 

X 


Turbulence  structural  coefficients  (equations  12-15 
Skin  friction  coefficient 
Specific  heat  at  constant  pressure 
Acceleration  parameter  (equation  5) 

Shape  factor,  5  /0 

Dissipation  length  scale  (equation  22) 

Turbulence  kinetic  energy  (equation  lb) 

Reynolds  number 

Reynolds  number  based  on  momentum  thickness 
Time  from  beginning  of  turbulent  period 
Turbulence  intensity 
Wall  temperature 
Freestream  temperature 

Instantaneous  value  of  velocity  components 

Fluctuating  velocity  components 

Mean  streamwise  velocity 

Dimensionless  velocity,  u/UT 

Friction  velocity 

Ensemble  averaged  velocity 

Velocity  fluctuations  (squared)  per  Hertz 

Distance  from  leading  edge 


Y 


Distance  from  wall 
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Subscripts 
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Dimensionless  distance,  YUT/v 
Intermit  tency 
Boundary  layer  thickness 
Displacement  thickness 

Long-term  average  boundary  layer  thickness 

Time  interval  from  beginning  to  end  of  a  turbulent  period 

Momentum  thickness 

Streamwise  integral  length  scale 

Kinematic  viscosity 


Freestream 
Fully  turbulent 
Upstream 
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Figure  3.  Oiagrams  of  the  Wedge  Inserts  Installed  in  the  Wind  Tunnel 
Test  Section 
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Figure  4.  Velocity  Distribution  Along  the  Flat  Test  Wall  with  the  2.3°  Angle 
Wedge  (Wedge  1)  and  No  Turbulence  Grid 
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Figure  6.  Turbulence  Intensity  and  Longitudinal  Integral  Scale  as  a  Function  of 
Speed  at  the  Entrance  to  the  Tunnel  Test  Section  —  Turbulence 

Grids  1,  2,  and  3 
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Figure  9.  Distributions  of  Total  Turbulence  Intensity  in  the  Test  Section  for 
K  =  0.20  x  10  -  6 
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Figure  11.  Distribution  of  the  Components  of  the  Turbulence  in  the  Test  Section 
for  K  =  0.75  x  10  ”  6  and  Grid  3 
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Figure  13.  Distributions  of  the  Streamwise  Integral  Length  Scale  in  the  Test 
Section  for  the  Four  Acceleration/Grid  Combinations 
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Figure  14.  U  &  V-Time  Records  Showing  a  Low  Freestream  Turbulence, 
Moderate  Intermittency  Flow,  Y/d  =  0.275 
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Figure  15.  Sample  U  &  W-Time  Records  Showing  a  Low  Freestream  Turbulence. 
Moderate  Intermittency  Flow,  Y Id  =  0.275 
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Figure  16.  Sample  U  &  V-Time  Records  Showing  a  Low  Freestream  Turbulence, 
High  Intermittency  Flow,  Y/d  =  0.183 
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Figure  17.  Sample  U  &  W-Time  Records  Showing  a  Low  Freestream  Turbulence, 
High  Intermittency  Flow,  Y/d  =  0.183 
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Figure  18.  Sample  UV,  UW  &  U-Time  Records  Showing  a  Low  Freestream 
Turbulence,  Low  Intermlttency  Flow,  Y/d  =  0.355 
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Figure  19.  Sample  UV,  UW,  &  U-Time  Records  Showing  a  Moderate  Freestream 
Turbulence,  Low  Intermittency  Flow,  Y Id  =  0.286 
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Figure  20.  Sample  UV,  UW,  &  U-TIme  Records  Showing  a  High  Freestream 
Turbulence,  Moderate  Intermittency  Flow,  Y/d  =  0.297 
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Figure  21.  Sample  UV,  UW,  &  U-Time  Records  Showing  a  Low  Freestream 
Turbulence,  High  Intermittency  Flow,  Y/d  =  0.145 
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Figure  22.  Sample  UV,  UW  &  U-Time  Records  Showing  a  Moderate  Freestream 
Turbulence,  High  Intermittency  Flow,  Y/d  =  0.172 
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Figure  23.  Sample  UV,  UW  &  U-Time  Records  Showing  a  High  Freestream 
Turbulence,  High  Intermittency  Flow,  Y/d  =  0.183 
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Figure  24.  Sample  U*Time  Records  at  3  Distances  From  the  Wall  for  a  Moderate 
Freestream  Turbulence,  Moderate  Intermittency  Flow 
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Figure  28C.  Profiles  of  Mean  and  Fluctuating  Quantities 
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Figure  28F.  Profiles  of  Mean  and  Fluctuating  Quantities 
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Figure  28G.  Profiles  of  Turbulent  Stresses 
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Figure  28K.  Profiles  of  Turbulent  Stresses 
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Figure  44.  Boundary  Layer  Spectral  Distributions 


Figure  45.  Boundary  Layer  Spectral  Distributions 
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Figure  53.  Predicted  Heat  Transfer-McDonald-Kreskovsky  Transition  Model 
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b)  STRONG  ACCELERATION  K  =  0  75E-06.  O  GRID  2.  A  GRID  3 

Figure  55.  Predicted  Skin  Friction-Cebeci-Smith  Model  with  Experimental  Intermittency 
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b)  STRONG  ACCELERATION  K  =  0  75E  06  O  NO  GRID  A  GRID  1.  +  GRID  2 

Figure  57.  Predicted  Heat  Transfer-McDonald-Kreskovsky  Model  with  Experimental 
Structural  Coefficients 
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-0.0820 
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• 
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29.4150 
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19.6100 

/V 

V 
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-0.0247 

-0.0740 

26.2730 
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-0.0123 

-0,0478 

30.7840 

39.0200 
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17.2091 
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-0.0024 
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32.2310 

57.6300 
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7.2037 
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Y/DELTA 

qll2/UEI»2 

qll2/UEII2 

qtl2/UEII2 

UV/UEU2 

UV/UEII2 

UV/UEII2 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

20NE 

0.2750 

0.008403 

0.005492 

0.014444 

-5.560E-04 

-7.096E-05 

-1.422E-03 

0.3250 

0.008131 

0.004950 

0.013944 

-5.339E-04 

-8.971E-05 

-1.343E-03 

0.3750 

0,007578 

0.003628 

0.013083 

-S.092E-04 

-2.480E-05 

-1.269E-03 

0.4250 

0.006724 

0.002857 

0.012031 

-4.564E-04 

-4.157E-05 

-1.09IE-03 

0.5000 

0.005591 

0.001715 

0.010627 

-4.266E-04 

-2.759E-05 

-9.662E-04 

0.6000 

0.004141 
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0.008799 

-3.072E-04 

-2.366E-05 

-7.50BE-04 

0.7000 

0.0023=5 

0.000467 

0.006438 

-1.835E-04 

-1.1176-05 

-4.976E-04 

0.8000 

0.001B28 

0.000354 

0.005182 

-1.273E-04 

-8. 177E-06 

-4.250E-04 

0.9000 

0.001310 

0.000272 

0.003812 

-8. 141E-05 

-3.320E-06 

-2.429E-04 

1.0000 

0.000925 

0.000272 

0.002937 

-6.647E-05 

-6.720E-06 

-2.173E-04 

1.1000 

0.000714 

0.000252 

0.002248 

-3.687E-05 

-3. 106E-06 

-1.075E-04 

5.0000 

0.000223 

0.000216 

0.000251 

-1.584E-07 

-5.090E-07 

5.939E-07 
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CROSS  HIRE  PROFILE  DATA 


1 

=  44.8 

6RID  NO. 

1 

K  *  0.2E-06 

Y/OELTA 

U7UE 

tr/UE 

U'/UE 

V7UE 

V/UE 

V/UE 

N/UE 

N/UE 

N/UE 

INTER- 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

0.1633 

0.0864 

0.0693 

0.0879 

0.0370 

0.0145 

0.0424 

0.0510 

0.0254 

0.0556 

C.2333 

0.0890 

0.0622 

0.0862 

0.0341 

0.0124 

0.0407 

0,0453 

0.0238 

0.0517 

0.2S33 

0.0892 

0.0582 

0.0867 

0.0334 

0.0117 

0.0386 

0.0432 

0.0211 

0.0486 

0.3333 

0.0863 

0.0455 

0,0857 

0.0314 

0.0101 

0.0368 

0.0390 

0.0171 

0.0466 

0.4167 

0.0765 

0.0356 

0.0782 

0.0280 

0.0083 

0.0336 

0.0349 

0.0142 

0.0429 

0.5000 

0.0652 

0.0261 

0.0685 

0.0248 

0.0074 

0.0302 

0.0330 

0.0137 

0.0386 

0.5833 

0.0544 
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0.0596 

0.0230 
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0.0282 

0.0114 

0.0347 
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0.0444 
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0.009! 

0.0251 

1.0000 
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0.0122 
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0.0091 
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1.1667 
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0.0176 

0.0095 

0.0068 

0.0113 
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0.0088 

0.0164 

3.3333 

0.0077 

0.0098 

0.0076 

0.0092 

0.0087 

0.0096 

0.0078 

0.0078 

0.0087 

T/OELTA  U'H2/qM2  U'»»2/q»l2  U'«2/q»»2  VII2/qtl2  V»t2/q»l2  V*I2/qII2  N'H2/qH2  N'»l2/qH2  N'H2/qH2 

INTER-  INTER-  INTER- 

COMPOSITE  TURBULENT  TURBULENT  COMPOSITE  TURBULENT  TURBULENT  COMPOSITE  TURBULENT  TURBUlENT 

ZONE  ZONE  Z  ONE  ZONE  ZONE  ZONE 

0.1833 

0.6536 

0.8497 

0.6127 

0.1200 

0.0374 

0.1430 

0.2263 

0.1129 

0.2443 

0.2333 

0.7117 

0.8457 

0.6429 

0.1051 

0.0330 

0.1373 

0.1832 

0.1212 

0.2198 

0.2833 

0.7274 

0.8551 

0.6613 

0.1023 

0.0341 

0.1319 

0.1703 

0.1106 

0.2069 

0.3333 

0.7483 

0.8403 

0.6761 

0.0995 

0.0415 

0.1246 

0.1522 

0.1182 

0.1990 

0.4167 

0.7451 

0.8247 

0.6733 

0.1003 

0.0453 

0.1250 

0.1546 

0.1300 

0.2017 

0.5000 

0.7141 

0.7402 

0.6618 

0.1035 

0.0585 

0.1292 

0.1823 

0.2013 

0.2091 

0.5833 

0.6919 

0.6817 

0.6446 

0.1236 

0.0865 

0.1372 

0.1845 

0.2316 

0.2180 

0.6667 

0.6621 

0.6557 

0.6272 

0.1205 

0.1183 

0.1342 

0.2174 

0.2261 

0.2386 

0.8333 

0.5949 

0.5654 

0.5794 

0.1679 

0.1737 

0.1509 

0.2372 

0.2608 

0.2697 

1.0000 

0.5222 

0.4579 

0.5281 

0.1928 

0.1934 

0.1529 

0.2850 

0.3486 

0.3190 

1.1667 

0.4132 

0.4688 

0.4453 

0.2522 

0.2000 

0.1786 

0.3346 

0.3312 

0.3761 

3.3333 

0.2905 

0.4140 

0.2578 

0.4J42 

0.3290 

0.4118 

0.2954 

0.2570 

0.3303 

Y/DELTA 

UV/Q**2 

UV/8442 

INTER- 

UV/8412 

6AMMA 

1 

6AHNA 

1 

COMPOSITE 

TURBULENT 

ZONE 

TURBULENT 

ZONE 

UV 

UV 

NV 

NV 

■scscssssxssssssssssssr: 

:::::;:znxuc;r 

0.1833 

-0.0941 

0.0037 

-0.1058 

71.8950 

24.2120 

66.1270 

19,4100 

0.2333 

-0.0863 

-0.0209 

-0.0986 

67.3080 

25.0130 

67.8710 

23.0120 

0.2833 

-0.0805 

-0.0224 

-0.0926 

71.9030 

23.8120 

74.6880 

21.6110 

0.3333 

-0.0800 

-0.0416 

-0.0877 

70.3380 

24.4130 

64.7980 

25.8130 

0.4167 

-0.0734 

-0.0433 

-0.0807 

67.1880 

27.8140 

61.0170 

26.2130 

0.5000 

-0.0646 

-0.0320 

-0.0748 

65.2820 

28.4150 

68.5250 

26.0130 

0.5833 

-0.0739 

-0.0446 

-0.0772 

67.7590 

27.0140 

59.5110 

27.0140 

0.6667 

-0.0652 

-0.0539 

-0.0699 

61.6140 

34.6180 

60.4070 

27.2140 

0.8333 

-0.0686 

-0.0339 

-0.0427 

68.0610 

38.0190 

49.6110 

32.8170 

1.0000 

-0.0571 

0.0034 

-0.0505 

63.4120 

46.0240 

38.6010 

35.0180 

1.1667 

-0.0348 

0.0041 

-0.0339 

54.6620 

54.2280 

22.1540 

29.4150 

3.3333 

0.0191 

0.0167 

0.0141 

45.7070 

64.0330 

2.4283 

9,0046 

Y/DELTA 

qll2/U£lt2 

ql42/UEI«2 

qll2/UEH2 

UV/UEM2 

UV/UEII2 

UV/UEII2 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

gxtsstszxcrssssttcxxxss: 

sxx&:xssx:sEss:iMtiuit«xt:sxs»ztKiuntmxs«stE»isixx: 

0.1833 

0.011427 

0.005672 

0.012621 

-1.075E-03 

2.085E-05 

-1.335E-03 

0,2333 

0.011138 

0.004660 

0.012091 

-9.618E-04 

-9.760E-05 

-1.193E-03 

0.2833 

0.010936 

0.004000 

0.011358 

-8.801E-04 

-8.951E-05 

-1.052E-03 

0.3333 

0.009947 

0.002459 

0.010870 

-7.957E-04 

-1.022E-04 

-9.536E-04 

0.4167 

0.007851 

0.001541 

0.009090 

-5.763E-04 

-6.679E-05 

-7.340E-04 

0.5000 

0.005958 

0.000928 

0.007097 

-3.846E-04 

-2.966E-05 

-5.307E-04 

0.5833 

0.004286 

0.000557 

0.005515 

-3.169E-04 

-2.4B3E-05 

-4.259E-04 

0,6667 

0.002991 

0.000531 

0.004112 

-1.949E-04 

-2.B59E-05 

-2. 6758-04 

0.8333 

0.001519 

0.000317 

0.002331 

-1.043E-04 

-1.076E-05 

-1 .462E-04 

1.0000 

0.000780 

0.000234 

0.001363 

-4.455E-05 

B.059E-07 

-6.877E-05 

1 . 1667 

0.000362 

0.000232 

0.000714 

-1.261E-05 

9.498E-07 

-2.423E-05 

3.3333 

0.000205 

0.000233 

0,000227 

3.915E-06 

3.900E-06 

3.I96E-06 
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TABLE  A8-A 


I  •  52.8 


CROSS  HIRE  PROPiLE  DATA 

GRID  NO.  i 


r  •  0.2E-06 


V/DELTA 

U'/UE  U'  /LIE 

U'/UE 

V/UE 

V'/UE 

V'/UE 

N'/UE 

N'/UE 

N'/UE 

INTER- 

INTER- 

INTER- 

COMPOSITE  TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

2  One 

ZONE 

ZONE 
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RO FILES  OF  MEAN  AND  ELUCTUATING  QUANTITIES 
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0.7682 
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0.6636 
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0.5476 
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0.0959 
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0.3565 

0.4292 
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0.3824 

0.4709 
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0.4247 

0.1881 
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0.2127 

0.4731 

0.4666 
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1.7143 

0.3101 

0.3343 

0.362: 

0.2278 

0.2166 

0.3806 

0.4621 

0,4419 

0.236° 
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0.3133 
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■ 
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0,325  : 
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C  .'.55 
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0.2-:  :•: 

0.5348 

0.54*3 
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0.06:: 

C .  r  : 

.  : .  z  z <  ; 

• 

(.24:: 

0.6'S: 
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0.06.7 
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Figure  A20-B 
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CROSS  MIRE  PROFILE  DATA 


* 


ft  S 


I 


p 

l 

*  36.8 

6RID  NO 

2 

K  »  0.75E-06 

1 1  l§ 

y/delta 

U/UE 

U'/UE 

U'/UE 

V  /UE 

V  /UE 

V'/UE  M'/UE 

M'/UE 

M  /uE 

INTER- 

INTER- 

INTER- 

■  m 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBU.ENT 

TURBulENT  COMPOSITE 

TURBUlENT 

TURBulENT 

1  | 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

0.2037 

0.0768 

0.0538 

0.0756 

0.0319 

0.0116 

0.0381  0.0424 

0.0246 

0.04*1 

IQS  iiC 

0.2963 

0.0703 

0.0479 

0.0676 

0.0286 

0.0134 

0.0332  0.0382 

0.0216 

0.0442 

K)  wC 

0.3889 

0.0605 

0.0317 

0.0602 

0.0256 

0.0096 

0.0303  0.0339 

0.0193 

0.036* 

HQ 

0.4815 

0.0486 

0.0225 

0.0509 

0.0218 

0.0077 

0.0260  0.0295 

0.0162 

0.0343 

IB 

0.5741 

0.0402 

0.0189 

0.0439 

0.0193 

0.0073 

0.0233  0.0269 

0.0172 

0.0316 

H  « 

0.7593 

0.0269 

0.0152 

0.0314 

0.0155 

0.0081 

0.0188  0,0223 

0.0172 

0.0262 

|M 

0.9444 

0.0170 

0.0127 

0.0226 

0.0119 

0.0086 

0.0146  0.0181 

0.0156 

0.0222 

Kq  ks 

1.1296 

0.0123 

0.0117 

0.0173 

0.0104 

0.0089 

0.0131  0.0156 

0.0146 

0.0196 

Is 

1.3148 

0.0104 

0.0112 

0.0148 

0.0101 

0.0097 

0.0127  0.0151 

0.0147 

0.0201 

Y/DELTA 

U’lt2/qtt2 

U’M2/qll2 

U'M2/q!!2 

V'tt2/qlI2 

VII2/q»l2  V'tt2/qt>2  M'U2/qU2  M  t!2/qt«2 

M  H2/qtl2 

B  cy 

INTER- 

INTER- 

INTER- 

f  4 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBUlENT  composite 

turbulent 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

is 

0.2037 

0.6768 

0.7940 

0.5967 

0.1161 

0,0363 

0.1504  0.2071 

0.16*6 

0.2530 

0.29o : 

0.6843 

0.7792 

0,6009 

0.1122 

0.0611 

0.142’  0.2033 

0.1557 

0.2562 

► 

0.3839 

0.6696 

0,6830 

0.5979 

0.1190 

0.0621 

0.1505  0.2115 

0.2545 

0.2516 

s 

0.4815 

0.6363 

0.5634 

0,5822 

0.1276 

0.0653 

0.1517  Q.23ei 

0.3713 

0.2661 

ff  v 

0.5741 

0.59e2 

0.5069 

0.5553 

0.1366 

0.0752 

0.1539  0.2672 

0.4180 

0.2906 

d  § 

0.7593 

0,4952 

0.3884 

0.4879 

0.1632 

0.1098 

0.1724  0.341c 

0.5016 

0.33s* 

0.9444 

0.3807 

0.3342 

0.4193 

0.1863 

0.1521 

0.1738  0.4330 

0.513' 

0.4069 

XI 

1.1296 

0.2988 

0.3170 

0,7V" 

0.2124 

0.1812 

0,1957  0.4898 

0.511- 

0.15': 

1  i 

1.3148 

0.2460 

0.2900 

0.257/ 

0.2327 

0.2124 

0.2003  0.5213 

0.457T 

0.512l 

R  * 

y/delta 

uv/qu: 

uv/on: 

Uv'/QHZ 

6AMMR 

4  GAMNfi 

f 

*/ 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

Uv 

UV  MV 

MV 

F£ 

ZONE 

ZONE 

It 

0.2037 

-0.0809 

-0.0252 

-0.0956 

66.2990 

24.4130  65.0900 

22.4110 

0.29s3 

-0.0843 

-0.0542 

-0.0977 

70.4330 

24.8130  66.3650 

22.6120 

IN 

G.3SS* 

-0.064c 

-0.0636 

-0,0984 

68.5300 

22.4110  66.6090 

23.812, 

Eh 

0.481: 

-0.0824 

-0.03-6' 

-0.0954 

67.2670 

25.0130  63.5860 

25.8130 

I Sr  V> 

0.5741 

-0.0762 

-0.0145 

-0.0857 

65.0200 

30.4160  57.8410 

26.81': 

Ec  'Tv 

0.7593 

-0.0779 

-0.00' 5 

-0.0887 

60.5070 

30.6160  50.8200 

30.8160 

0.9444 

-0.0514 

-0.0025 

-0.0662 

47.9640 

42.8220  32.1570 

37.4190 

P 

1.1296 

-0.0232 

0.0057 

-0.0496 

30.7990 

45.0230  16.2140 

33.8170 

RS  \i 

1.3148 

0.0041 

0.0109 

-0,0135 

19.5770 

49.8260  5.8940 

17.2090 

Y/DELTA 

qll2/UEII2 

qtl2/UEM2 

q«*2/UE*  *2 

UV/UEH2 

UV/UEII2 

UV/UEII2 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

0.2037 

0.008718 

0.003662 

0.009592 

-7.051E-04 

•SSSSSCSSSSS 

-9.245E-05 

-9.193E-04 

0.29e3 

0.007240 

0.002940 

0.007669 

-6.102E-04 

-1.592E-04 

-7.492E-04 

0.386* 

0.0054TZ 

0.001474 

0.00605* 

-4.630E-04 

-9.379E-05 

-5.964E-04 

0.4815 

0,003715 

0,000899 

0.004445 

-3.061E-04 

-3.296E-05 

-4.240E-04 

0.5741 

0.002722 

0.000714 

0.0034*0 

-2.074E-04 

-1.064E-05 

-2.991E-04 

0.75*3 

0.00145* 

0.0005*5 

0.00202* 

-1.137E-04 

-5.652E-06 

-1.799E-01 

0.9145 

0.00075* 

0.0004*9 

0.001222 

-3.900E-05 

-1.1B9E-06 

-B.089E-05 

1.12*6 

0.0005-.: 

C. 000430 

0. 000865 

-1.168E-05 

2.468E-06 

-4.292E-C5 

1.3149 

0.000439 

0.000436 

0.000797 

1.802E-06 

4.756E-06 

-1.079E-05 
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Figure  A21-B 


X  *  48. e 


SINGLE  HIRE  PROFILE  DATA 


GRID  NO.  2 


K  *  0.75E-06 


U/uE  U/UE  U/UE  IT/UE  U  /HE  U'/UE 


COMPOSITE 

TURBUlENT 

ZONE 

TURBULENT 

ZONE 

COMPOSITE 

TURBULENT 

ZONE 

TURBULENT 

ZONE 

u 

u 

0.0469 

0.7351 

0.7162 

0.7359 

0.0753 

0.0729 

0.0753 

95.8120 

7.4038 

0.0500 

0.7361 

0.7021 

0,7383 

0.0745 

0.0601 

0.0748 

94.2060 

6.0031 

0.0531 

0.7471 

0.7423 

0.7473 

0.0717 

0.0583 

0.0722 

96.0480 

5.002s 

0.0563 

0.7511 

0,7533 

0.7510 

0.0707 

0.0586 

0.0715 

93.4020 

7. 603* 

0.0594 

0.7598 

0.7721 

0.7589 

0.0699 

0.0554 

0.0707 

93.4140 

8.4043 

0.0625 

0.7650 

0.7872 

0.7639 

0,0661 

0.0557 

0.0685 

95.0080 

7.4036 

0.0700 

0.7720 

0.8090 

0.7695 

0.0676 

0.0698 

0.0669 

94.3010 

6.4033 

0.0750 

0.7779 

0.8128 

0.7761 

0.0659 

0.0655 

0.0654 

95.1920 

7.4038 

0.0813 

0.7913 

0.8357 

0.7871 

0.0656 

0.0647 

0.0643 

92.7590 

8,6050 

0.0938 

0.8039 

O.S<7C 

0.8002 

0.0640 

0.0644 

0.062? 

94.02''.' 

6.0 

0.1094 

0.8194 

0.8602 

0.8162 

0.0616 

0.0738 

0.0598 

96  .OebC 

7 . 2 «.  7  ~ 

0.1250 

0.8290 

0 . 8  fc  S  C1 

0.8256 

0.0597 

0.0733 

0.0575 

96.7210 

5.012c 

0.140s 

0.8469 

0.89(1 

0.8397 

0.0613 

0.0718 

0.0565 

93.089/ 

7.037c 

0. 15t3 

0.8568 

0.8211 

0.8525 

0.05«9 

0.0680 

0.0554 

93.3220 

6.4,4'' 

0.1875 

0.6751 

0,93s4 

0.8663 

0.0579 

0.0650 

0.0535 

Vi . 9fc6 v 

6.604* 

0.23<4 

0.897c 

0.9;47 

0.892s 

0.0517 

0.0481 

0.0465 

94.0c31 

6  .  C  0  *4 

0.3125 

0.9229 

0.9730 

0.9184 

0.044? 

0.0396 

0.0427 

92.75*0 

9,004c 

0.390c 

0.94;; 

0.9’53 

0.843; 

0.03«! 

0.0353 

0.0*76 

87.2770 

1 3 ,  c ;  -  / 

0.4688 

0.9619 

0.9804 

0.9564 

0.0335 

0.0410 

0.0330 

84.429': 

21.8111 

0.6250 

0.9810 

0.8803 

0.8749 

0.0220 

0.0155 

0.0247 

58.6030 

62.071: 

0.7613 

0.9911 

o,cc:u 

0.9831 

0.0145 

0.0124 

0.0196 

24.2260' 

55.6283 

0.9375 

0,99;: 

0.996s 

0.9661 

0.009' 

0.0102 

0.017c 

5.308  = 

i; .  ;  : 

1.093s 

0.9963 

0.9853 

0.9935 

0.0071 

0.008c 

0.0142 

1.31®2 

*.0C3t 

1.2500 

i.oo 

1.0000 

1.005.7 

0 . Out  1 

0.0080 

O.OlOE 

0.1511 

0.60:7 

1.4063 

1.0011 

l  ,0011 

0,98'; 

0.0062 

o.ooec 

0.0048 

0.06*5 

O.e  y.  3 

2.5000 

l.ooo : 

1.0003 

0.9846 

0.005s 

0.0056 

0.0000 

0.0000 

0.0000 

TABLE  A21-A 
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CROSS  HIRE  PROFILE  DATA 


l  «  48.8 


GRID  MO.  2 


K  =  0.75E-06 


Y/DELTA 

U7UE 

U7UE 

U7UE 

V7UE 

V7UE 

V7UE 

N7UE 

N/UE 

N/UE 

INTER- 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

Z  ONE 

ZONE 

ZONE 

0.1719 

0.0655 

0.0723 

0.0619 

0.0315 

0.0177 

0.0322 

0.0413 

0.0312 

0.0424 

0.2500 

0.0578 

0.0570 

0.0554 

0.0279 

0.0126 

0.0286 

0,0365 

0.0300 

0.03T3 

0.3281 

0.0496 

0.0433 

0.0484 

0.0245 

0.0085 

0.0252 

0.0330 

0.0311 

0.0333 

0.4063 

0.0423 

0.0425 

0.0416 

0.0216 

0.0180 

0.0221 

0.0284 

0.0220 

0.0289 

0.4644 

0.0352 

0.0332 

0.0353 

0.0186 

0.0122 

0.0191 

0.0249 

0.0199 

0.0255 

0.6406 

0.0241 

0.0221 

0.0246 

0.0140 

0.0091 

0.0143 

0.0199 

0.0150 

0.0205 

0.7969 

0.0157 

0.0184 

0.0163 

0.0106 

0.0070 

0.0107 

0.0160 

0.0120 

0.0167 

0.9531 

0.0108 

0.0141 

0.0118 

0.0090 

0.0076 

0.0091 

0.0135 

0.0114 

0.0147 

1.1094 

0.0083 

0.0107 

0.0093 

0.0063 

0.0077 

0.0084 

0.0121 

0.0114 

0.0134 

1.2656 

0.0071 

0.0102 

0.0079 

0.0084 

0.0075 

0.0087 

0.0114 

0.0111 

0.0124 

1.4219 

0.0068 

0,0098 

0,0071 

0.0087 

0.0080 

0.0089 

0.0108 

0.0105 

0.0128 

2.5156 

0.0066 

0.0101 

0.0088 

0.0108 

0.0098 

0.0091 

0,0097 

0.0095 

0.0111 

T/DELTA 

U'H2/qlt2 

UU2/qU2 

U'U2/qll2 

V’!t2/qtt2 

V»l2/ql»2 

V'!t2/qt!2 

N'H2/q»l2 

rtt2/qtl2 

N'H2/qll2 

INTER- 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

tttiixstxtst 

:xx:s>xsusu«scxisc3xss:3xs:sxtxicss:s::s:3ss8: 

::zxx:t:c::: 

sxstsscsxsss 

0.1719 

0.6140 

0.8027 

0.5754 

0.1400 

0.0472 

0.1535 

0.2459 

0.1501 

0.2711 

0.2500 

0.6127 

0.7558 

0.5812 

0.1410 

0.0357 

0.1528 

0.2463 

0.2085 

0.2660 

0.3281 

0.5928 

0.6609 

0.5725 

0.1430 

0.0228 

0.1534 

0.2642 

0.3163 

0.2741 

0.4063 

0.5844 

0.6911 

0.5688 

0.1504 

0.1217 

0.1566 

0.2652 

0.1872 

0.2745 

0.4644 

0.5619 

0.6678 

0.5512 

0.1547 

0.0894 

0.1583 

0.2834 

0.2428 

0.2905 

0.6406 

0.4961 

0.6106 

0.4934 

0.1640 

0.1028 

0.1635 

0.3399 

0.2866 

0.3431 

0.7969 

0.3999 

0.6436 

0.4039 

0.1786 

0.0893 

0.1713 

0.4215 

0.2670 

0.424B 

0.9531 

0.3076 

0.5236 

0.3178 

0.2102 

8.1454 

0.1669 

0.4822 

0.3311 

0.4953 

1.1094 

0.2401 

0.4006 

0.2585 

0.2378 

0.1860 

0.2067 

0.5221 

0.4134 

0.5348 

1.2656 

0.1990 

0.3700 

0.2127 

0.2784 

0.1953 

0.2555 

0.5226 

0.4347 

0.5318 

I. 4219 

0.1957 

0.3591 

0.1725 

0.3113 

0.2318 

0.2682 

0.4930 

0.4090 

0.5593 

2.5156 

0.1715 

0.3572 

0.2834 

0.4526 

0.3282 

0.2866 

0.375B 

0.3146 

0.4301 

Y/DELTA 

UV/BU2 

UV/8II2 

INTER- 

UV/BII2 

6AMMA 

< 

6AMMA 

f 

COMPOSITE 

TURBULENT 

ZONE 

TURBULENT 

ZONE 

UV 

UV 

NV 

NV 

xxsxxxxxcxxxixxxxxxxxxxxxxxxxsKxxxxxxxxtxrxxxstxistcxxxsxtxxxeszsxcxxxs:: 

SSSSSEZSSSSSS 

xrr:s:::x: 

0.1719 

-0.1084 

-0.0524 

-0.1158 

94.8870 

6.2032 

93.1510 

7.4038 

0.2500 

-0.1075 

-0.0445 

-0.1147 

94.7690 

9.2047 

93.6630 

7.0036 

0.3281 

-0.1100 

-0.0077 

-0.1174 

93.9930 

8.6044 

95.8710 

3.8019 

0.4063 

-0.1065 

-0.0682 

-0.1132 

94.8920 

8.2042 

93.1860 

10.0050 

0.4844 

-0.1039 

-0.0612 

-0.1069 

94.2670 

10.6050 

91.7670 

11.0060 

0.6406 

-0.0933 

-0.0241 

-0.0951 

94.1570 

12.2060 

88.3630 

12.4060 

0.7969 

-0.0717 

-0.0130 

-0.0709 

94.0700 

13.4070 

82.8280 

27.2140 

0.9531 

-0.0411 

0.0041 

-0.0413 

89.5540 

26.6140 

59.4540 

51.4260 

1.1094 

-0.0100 

0.0084 

-0.0135 

82.7020 

43.8220 

34.8280 

55.8290 

1.2656 

0.0156 

0.0104 

0.0122 

72.3030 

59.2300 

20.1050 

40.6210 

1.4219 

0.0189 

0.0163 

0.0124 

67.8560 

59,4300 

11.6090 

31.2160 

2.5156 

0.0324 

0.0120 

-0.0633 

73.2040 

59.2300 

10.4710 

26.6140 

f 


l 

i 

i 


■. 

a 

« 

<* 


Y/DELTA 

qtl2/UEH2 

qlt2/U£H2 

qtt2/UEH2 

UV/UEM2 

UV/UEII2 

UV/UEM2 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

0.1719 

0.007005 

0.006549 

0.006695 

-7.591E-04 

-3.429E-04 

-7.752E-04 

0.2500 

0.005461 

0.004379 

0.005283 

-5.870E-O4 

-1.949E-04 

-6.061E-04 

0.3281 

0.004161 

0.003101 

0.004103 

-4.576E-04 

-2.391E-05 

-4.817E-04 

0.4063 

0.003067 

0.002616 

0.003083 

-3.268E-04 

-1.784E-04 

-3.49IE-04 

0.4844 

0.002215 

0.001648 

0.002267 

-2.S01E-04 

-1.009E-04 

-2.422E-04 

0.6406 

0.001175 

0.000799 

0.001234 

-I.096E-04 

-1.92IE-0S 

-1.174E-04 

0.7969 

0.000618 

0.000547 

0.000667 

-4.432E-05 

-7.105E-06 

-4 . 730E-05 

0.9531 

0.000382 

0.000394 

0.000442 

-1.568E-05 

1.61BE-06 

-1.826E-05 

1.1094 

0.000286 

0.000318 

0.000337 

-2.87IE-06 

2.680E-06 

-4.544E-06 

1.2656 

0.000252 

0.000287 

0.000294 

3.917E-06 

2.989E-06 

3.577E-06 

1.4219 

0.000240 

0.000273 

0,000295 

4.519E-06 

4.429E-06 

3.652E-06 

2.5156 

0,000253 

0.000290 

0.000288 

B.179E-06 

3.487E-06 

-1.826E-05 

TABLE  A21-B 
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PROFILES  OF  MEAN  AND  FLUCTUATING  QUANTITIES 


Y/DELTA  Y/DELTA  Y/DELTA 

Figure  A22 


I  --  4.6 


K  =  0.75E-06 


SINGLE  HIRE  PROFILE  DATA 

GRID  NO.  3 


I/DElTA  y/UE  U/UE 


U/UE  U'/UE  U/UE 


U/UE  6ANNA 


COMPOSITE 

INTER- 

TURBULENT 

ZONE 

0.1678 

0.3892 

0.3797 

0.1709 

0.4177 

0.4084 

0.1933 

0.4446 

0.4357 

0.2000 

0.4602 

0.4506 

0.2089 

0.4027 

0.4742 

0.2333 

0.5162 

0.5079 

0.2344 

0.5285 

0.5215 

0.2456 

0.5505 

0.5421 

2.6556 

0.5912 

0.5826 

0.2933 

0.6396 

0.6329 

0.3089 

0.6603 

0.6538 

0.3311 

0.6920 

0.6899 

0. 3889 

0.7544 

0.7534 

0.4444 

0.8183 

0.8225 

0.5000 

0.8604 

0.864o 

0.8;23 

0.8975 

0.6111 

0.9130 

0.9214 

0,66e' 

0.9297 

0.9349 

0.7778 

0.9637 

0.9683 

0.8889 

0.9772 

0.9804 

1.1111 

0.9959 

0.9974 

1.3333 

0.9991 

0.999' 

1.5556 

1.0022 

1.0019 

2.7778 

0.9966 

0.9983 

INTER- 


TURBULENT 

ZONE 

CGNFOSITE 

TURBULENT 

ZONE 

0.5153 

0.0987 

0.0867 

0.5384 

0.1026 

0.0912 

0.5730 

0.1013 

0.0911 

0.5770 

0.106B 

0.0968 

0.5762 

0.1066 

0.0981 

0.5985 

0.1115 

0.1024 

0.5954 

0,1089 

0.1023 

0.6341 

0.1128 

0.1041 

0.6685 

0.1152 

0.1074 

0.6899 

0.1127 

0.1057 

0.7138 

0.1127 

0.1071 

0.7098 

0.1138 

0.1085 

0.7604 

0.1128 

0.1073 

0.7945 

0.1086 

0.1024 

0.8355 

0.1005 

0.0939 

0.8575 

0.0973 

0.090c 

0.8657 

0.0936 

0.0850 

0.8895 

0.0865 

0.0788 

0.9249 

0.0764 

0.0699 

0.9456 

0.0705 

0.0647 

0.9729 

0.0646 

0.0621 

O.9035 

0.0613 

0.0604 

1.0126 

0.0585 

0.0582 

1.0546 

0.0536 

0.0537 

TURBULENT 

U 

Li 

ZONE 

0.1485 

6.5900 

16.8090 

0.1534 

6.6598 

17.6090 

0.1449 

6.0605 

16.4080 

0.1467 

7.2157 

19,0100 

0.1445 

7.3356 

20.2100 

0.1564 

8.7500 

24.6130 

0.1428 

8.9421 

23.8120 

0.1552 

8.5015 

23,0120 

0.1492 

9.4826 

25.2130 

0.1461 

11.0710 

29.4150 

0.1396 

10.34B0 

26. e 150 

0.1502 

10.0790 

26,61*0 

0.1422 

13.2680 

37, 41  =  .. 

0.1361 

14.21U1 

39,4200 

0.1302 

13.1070 

36.415., 

0.1291 

12.2260 

3c. 0180 

0.1224 

14.1110 

36.2200 

0.1243 

10.9070 

35.4180 

0.1119 

9.6030 

32 . 01 eO 

0.1085 

8.5528 

26.015; 

0.0936 

5.7454 

19.6100 

0.0814 

3. 7141 

13,607.; 

0.0741 

2.2490 

8.60-45 

0.0549 

0.26«0 

1.2006 

TABLE  A22 
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COMPOSITE 


PROFILES  OF  TURBULENT  STRESSES 

X  =  8.8  K  =  0.75  E— 6  TE  =  4.4  * 


Y/DELTA  Y/DELTA  Y/DELTA 

Figure  A23-B 


I, 


SINSlE  HIRE  PROFILE  DATA 

X 

=  6.8 

6RID  NO. 

3 

K  «  0.75E-06 

U/UE 

U'UE 

U/UE 

U/UE 

U/UE 

U/UE 

COMPOSITE 

INTEP- 

TURB'JLENT 

TURBULENT 

COMPOSITE 

INTER- 

TURBUlENT 

TURBULENT 

TABLE  A23-A 


W 


CROSS  HIRE  PROFILE  DATA 


l  ■-  8.8 


GRIB  NO.  3 


K  *  0.75E-06 


Y/  DELTA 


L1  /'JE  U  /OE  V  /UE  V'/UE 


k  /UE  k  /UE 


INTER-  INTER-  INTER  - 

COMPOSITE  TURBUlEN7  TURBU-EN'  COMPOSITE  TURBUlENt  TURBUiENt  C0HP05ITE  TURBUlEN7  TuRBu-EN' 

ZONE  ZONE  ZONE  ZONE  ZONE  ZONE 


Y/DELTA  U'112/qttl 

U  U2/q»l2 

U  t *2 / q 1  * 2 

V  **2/q» 12 

VII2/qH2 

V'»l2/q«»2 

k'll2/qM2 

N'lt2/qtl2 

k'U2/q»l2 

COMPOSITE 

INTER- 

TURBULENT 

TURBULEN7 

COMPOSITE 

INTER- 

TURBULENT 

TURBULENT 

COMPOSITE 

INTER- 

TURBULENT 

TURBuiEn’ 

COMPOSITE 

TURBULENT 

ZONE 

TURBULENT 

ZONE 

uv 

UV 

kv 

kv 

0.2972 

-0.0654 

-0.0402 

-0.0858 

42.2260 

49.0250 

39.7260 

4t .  62* ; 

0.3789 

-0.062? 

-0.0*69 

-0.0767 

41.3830 

48.4250 

40.3640 

47.C24. 

0.459! 

-0.0695 

-0.0442 

-0,0878 

40.1540 

49.6250 

40.3590 

49.4250 

0.54c! 

-0.0528 

-0.03(5 

-0.0703 

42.1210 

56.6290 

40.0540 

52.02*0 

0.6216 

-0.0464 

-0.0226 

-0.0626 

42.6490 

51.8270 

38.8550 

48.0250' 

0.702’ 

-0.0440 

-0.0126 

-0.0604 

46.5470 

61.0310 

32.5850 

49.2250 

0.7838 

-0.0348 

-0.0071 

-0.0527 

42.2690 

61.2310 

32.4180 

47.8240 

0.864? 

-0.0254 

-0.0085 

-0.0386 

43.9340 

56.6300 

29.5470 

47.6240 

1,0811 

-0.0112 

0.0050 

-0.0244 

41.6240 

61.0310 

21.8880 

40.6210 

2.4324 

0.0043 

0.0058 

0.0010 

37.5820 

64.8330 

9.8950 

26.8140 

Y/DELTA  q»I2/UEU2  q»«2/UEII2  qtl2/UEl»2  UV/UEII2  UV/UEt 1 2  UV/UEII2 

INTER-  INTER- 

COMPOSITE  TURBULENT  TURBULENT  COMPOSITE  TURBULENT  TURBULENT 

ZONE  ZONE  ZONE  ZONE 

::si::::::::::::rixrz:ztxtrsstzt2Sxriixttiiit:tscxz:ziri:cxtttitii£szKsrr:::::::: 


0.2»73 

0.016900 

0.014712 

0.020333 

-1.106E-03 

-5.921E-04 

-1.745E-03 

0.378* 

0.014993 

0.011785 

0.018818 

-9.586E-04 

-5.529E-04 

-1.444E-03 

0.4595 

0.012931 

0.009559 

0.017176 

-B.967E-04 

-4.238E-04 

-1.5Q7E-03 

0.5405 

0.011884 

0.009744 

0.015423 

-6 . 389E-04 

-2 . 972E-04 

-1.084E-C 

0.6216 

O.OtlCl* 

0.008*96 

0.014392 

-5. 113E-04 

-1.920E-04 

-9.005E-0* 

0.702’ 

0.0094; i 

0.007272 

0.012706 

-4, 140E-04 

-9.136E-05 

-7.660E-0* 

0.7838 

0.009217 

0.007224 

0.012649 

-3.204E-04 

-5.108E-05 

-6 . 660E-0* 

O.Ssl® 

0.00890? 

0. 0075t5 

0.011378 

-2.261E-04 

-6.460E-05 

-4.387E-0* 

1.0811 

0.007886 

0.006823 

0.010577 

-B.829E-05 

3.413E-05 

-2.576E-04 

2.432* 

0.0068?2 

0.006420 

0.008715 

2.944E-05 

3 . 700E-05 

8.800E-06 

TABLE  A23-B 
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PROFILES  OF  MEAN  AND  FLUCTUATING  QUANTITIES 


Y/DELTA  Y/DELTA  Y/DELTA 

Figure  A24-B 


SIN6LE  WIRE  PROFILE  DATA 


l  *  12. B 


6RIC  WO.  3 


K  *  0.75E-06 


TABLE  A24-A 


V/DELTA 

U'UE 

COMPOSITE 

U/UE 

INTEA- 

TURBULENT 

ZONE 

U/UE 

TURBULENT 

ZONE 

U’/UE 

COMPOSITE 

U  /UE 

INTER- 

TURBULENT 

ZONE 

U'/UE 

TURBULENT 

ZONE 

6ANNA 

U 

1 

U 

3  3 

) 

A 

„N 

0.0600 

0.5003 

0.40*0 

0.5377 

0.1382 

0.1077 

0.1305 

70.9120 

*6.22*0 

^  r< 

0.0632 

0.5310 

0.4*36 

0.5580 

0.1353 

0.1109 

0.1306 

75.4790 

*7.02*0 

0.0668 

0.5376 

0.4535 

0.5698 

0,1363 

0.1086 

0.1320 

71.4930 

*7.02*0 

7  3 

0.0708 

0.5619 

0.4911 

0.5859 

0.13*4 

0.1137 

0.132* 

73.8110 

*4.2230 

V 

0.07*0 

0.5769 

0.5005 

0.60*4 

0.13*8 

0.11*5 

0.1J09 

72.9100 

*2.0220 

0.0772 

0.6002 

0.5258 

0.6202 

0.1322 

0.1137 

0.1297 

78.2*30 

38.6200 

0.0820 

0.6108 

0.55*0 

0.6329 

0.1315 

0.1191 

0.1295 

71.3270 

*2.6220 

S-  j 

0.0860 

0.6388 

0.595* 

0.651* 

0.1271 

0.1162 

0.127* 

76.7210 

35.2180 

0.09*0 

0.6630 

0.6202 

0.67*9 

0.1267 

0.1182 

0.126* 

77.4720 

39.0200 

0.1020 

0.6815 

0.6391 

0.6927 

0.1237 

0.1220 

0.1217 

78.3020 

39.4200 

$ 

0.1084 

0.70*1 

0.6827 

0.7103 

0.1196 

0.11*4 

0.120* 

77.2620 

34.8180 

A! 

0.1176 

0.7199 

0.7065 

0.7237 

0.1192 

0.1139 

0.120* 

76.9310 

36.0180 

0.1360 

0.7*88 

0.7499 

0.7*85 

0.1089 

0.103* 

0.1103 

78.0050 

38.2200 

y.  & 

0.1556 

0.7782 

0.7973 

0.7729 

0.1072 

0.1059 

0.1070 

77.7*60 

35.2180 

jtv.  ►! 

j 

0.1760 

0.7999 

0.8271 

0.7926 

0.10*6 

0.0987 

0.10*9 

77.9230 

39.2200 

A 

0.1952 

0.816* 

0.8*61 

0.8077 

0.1029 

0.0966 

0.1026 

76.3860 

*0.0200 

0.2156 

0.8327 

0.8611 

0.82*2 

0.0971 

0.1023 

0.096* 

77.7*80 

*0.0200 

•S  ' 

0.23*0 

0.8*59 

0.8826 

0.8353 

0.09*2 

0.0822 

0.09*7 

76.8080 

*0,0200 

i. 

0.2752 

0.8691 

0.90*5 

0.8560 

0.09vb 

0.0753 

0.0926 

72.8380 

*9.4250 

•j 

0.3152 

0.8873 

0.9286 

0.6720 

0.0867 

0.0710 

0.0872 

72.2030 

*5.6230 

■ "  " li 

1 

0.39*8 

0.9155 

0.9*26 

0.9038 

0.0799 

0.0683 

0.0818 

68.5960 

56.2300 

. .  ~ 

0.*9*< 

0.9*5’ 

0.962’ 

0,9326 

0.0711 

0.05E6 

0.0773 

5*.9E1C 

56.4300 

•j 

0.5932 

0.9628 

0.9726 

0.9517 

0.06*8 

0.0558 

0.072* 

*5.3360 

67,03*i 

r.  "J 

0.6936 

0.9800 

0.96*= 

0.9705 

0.058* 

0.0519 

0.0687 

33.38=0 

59.6310 

-i 

0.7932 

0.9855 

0.989’ 

0.98*5 

0.0532 

0.0*89 

0.0660 

22.0360 

*9.8260 

0.993s 

0.992’ 

0.9917 

0.9955 

0.0*85 

0.0540 

0.059* 

10.026: 

29.6151 

1.1932 

0.9®E6 

0.9961 

1.013* 

O.O**7 

0.0*39 

0.0607 

4.2361 

15.4080 

2.192* 

1.004* 

1.003* 

0.983* 

0.0*16 

0.0521 

0.0228 

0.164* 

1.0005 

■  •: 
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CROSS  MIRE  PROf JLE  ORTA 

GRID  NO.  3 


r  *  0.75C-06 


T/OElTA 

U'/UE 

U'/UE 

U'/UE 

V'/UE 

V'/UE 

V'/UE 

N'/UE 

N'/UE 

N'/UE 

INTER- 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

20NE 

0.2200 

0.0932 

0.0947 

0.0925 

0.0433 

0.0251 

0.0472 

0.0762 

0.0656 

0.0787 

0.2000 

0.0874 

0.0814 

0.0874 

0.0407 

0.0241 

0.0438 

0.0717 

0.0618 

0.0743 

0.4000 

0.0740 

0.0411 

0.0782 

0.0344 

0.0209 

0.0394 

0.0698 

0.0606 

0.0723 

0.5200 

0.0492 

0.0540 

0.0727 

0.0332 

0.0205 

0.0342 

0.0701 

0.0640 

0.0724 

0.5800 

0.0440 

0.0537 

0.0700 

0.0322 

0.0205 

0.0354 

0.0651 

0.0594 

0.0675 

0.6400 

0.0414 

0.0513 

0.0459 

0.0308 

0.0212 

0.0341 

0.0638 

0.0556 

0.0686 

0.7000 

0.0584 

0.0497 

0.0433 

0.0298 

0.0215 

0.0331 

0.0653 

0.0598 

0.0692 

0.0000 

0.0553 

0.0485 

0.0404 

0.0263 

0.0223 

0.0313 

0.0634 

0.0561 

0.0703 

0.4000 

0.0531 

0.0499 

0.0588 

0.0285 

0.0230 

0.0316 

0.0419 

0.0585 

0.0666 

1.0000 

0.0491 

0.0452 

0.0541 

0.0281 

0.0232 

0.0315 

0.0618 

0.0577 

0.0688 

1.2000 

0.0475 

0.0459 

0.0504 

0.0291 

0.0258 

0.0324 

0.0557 

0.0530 

0.0629 

2.2000 

0.0431 

0.0473 

0.0444 

0.0368 

0.0338 

0.0396 

0.0526 

0.0520 

0.0572 

r/DELTfl 

U’tl2/qU2 

Ut»2/qt»2 

U'U2/qt(2 

VU2/qH2 

V'tt2/Qlt2 

V»2/qll2 

N‘M2/qtl2 

Nt»2/ql«2 

Ntl2/qtt2 

INTER- 

INTER- 

INTER- 

COMPOS  I TE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

iss£tss:st:sssttutittttstet:tmss: 

ssssetctcrs: 

CSCtSCttttltlttSSItBCtSS 

0.2200 

0.5286 

0.6487 

0.5015 

0.1108 

0.0430 

0.1273 

0.3606 

0.3083 

0.3712 

0.2800 

0.5274 

0.5986 

0.5048 

0.1109 

0.0508 

0.1231 

0.3617 

0.3506 

0.3722 

0.4000 

0.4795 

0.4746 

0.4712 

0.1071 

0.0536 

0.1167 

0.4134 

0.4719 

0.4121 

0.5200 

0.4403 

0.4064 

0.4430 

0.0985 

0.0527 

0.1072 

0.4612 

0.5409 

0.4499 

0.5800 

0.4494 

0.4179 

0.4550 

0.1043 

0.0595 

0.1130 

0.4463 

0.5226 

0.4319 

0.6400 

0.4249 

0.4223 

0.4227 

0.1045 

0.0705 

0.1098 

0.4705 

0.5072 

0.4675 

0.7000 

0.3945 

0.3751 

0.4019 

0.1000 

0.0686 

0.1069 

0.5054 

0.5563 

0.4911 

0.8000 

0.3844 

0.3886 

0.3777 

0.0980 

0.0797 

0.0990 

0.5176 

0.5316 

0.5232 

0.9000 

0.3740 

0.3853 

0.3860 

0.1054 

©  .0789 

0.1084 

0.5206 

0.5358 

0.5056 

1.0000 

0.3408 

0.3413 

0.3367 

0.1088 

0.0880 

0.1107 

0.5504 

0.5707 

0.5526 

1.2000 

0.3612 

0.3744 

0.3349 

0.1322 

0.1150 
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IjNE 

;  >: 

0.1279 

0.07tl 

0.0772 

0.0.763 

0.042c 

C.0422 

0.0428 

0.0625 

O.OcOl 

-  - 

'V 

o.i etc 

0.0905 

0.0749 

0.0703 

0.0394 

0.0393 

0.0394 

0.056t 

C.0547 

0,(569 

0.2*41 

0.0666 

O.Ooe' 

0.0669 

0.0365 

0.03c2 

0.0366 

0.0602 

O.Oel" 

0,0596 

V  . 

0.3023 

0.0623 

0.0640 

0.0c27 

0.0341 

0.0341 

0.0341 

0.0563 

0.0546 

0.059. 

0.3605 

0.05S*> 

0.0573 

0.0586 

0.0319 

0.0292 

0.0320 

0.0544 

0.049c 

0.054' 

■ 

0,4767 

0.0514 

0.0419 

0.0524 

0.0279 

0.0227 

0.02B2 

0.0507 

0.0426 

0,051' 

0.5930 

0.0462 

0.0380 

0,0481 

0.0259 

0.0195 

0.0264 

0.0505 

0.0460 

0.0516 

0.6256 

0.0391 

0.0329 

0.0401 

0.0236 

0.0183 

0.0247 

0.0480 

0.0457 

0.053 ; 

0,9419 

0.C357 

0.0321 

0.0387 

0.0249 

0.0146 

0.0262 

0.0444 

0.040e 

0.04=4 

r  - 

fcV 

2.325c 

0.0282 

0.0267 

0.0286 

0.0334 

0.0000 

0.0346 

0.0352 

0,034'.' 

0.03=3 

.  / 

*  . 

Y/DElTA 

UII2/QM2 

U  H2/q»»2 

U  *  1 2  /  q  1 1 2 

V  t»2/qt!2 

V  1 1 2' q  *  1 2 

V  UZ  qli: 

N  111  ql«: 

M  112/qli: 

N  112/qlt: 

* 

INTER- 

IN’EA- 

INTER- 

r£ 

C0X90EITE 

TURBU.EN' 

TURBulEnt 

COMPOSITE 

TURBUlEn’ 

TURBULENT 

COMPOSITE 

TUR8UlEN7 

TURB--vENT 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

ZONE 

S 

0.1279 

0.5:25 

0.525* 

0.5010 

0.1563 

0.  IS*-” 

O.l5o4 

0.3412 

0.3194 

0.3426 

i 

o.  iSc : 

0.4969 

0.5574 

0.4957 

0.1542 

0.1514 

0.1542 

0.3409 

0.2=12 

0.2442 

C.4717 

0.459c 

0.4755 

0.1406 

0.1356 

0.1410 

0.3877 

0.404c 

0.3835 

0.3023 

0.456c 

0.497(, 

0.4542 

0.13c5 

0.139! 

0,1353 

0.4047 

0.3c3e 

0.415 

K 

0 . 3 1 . 5 

0.4t24 

G.49'7 

0.461 7 

0.1356 

0.127 7 

0.1355 

0.4016 

0.374c 

0.4026 

5b- 

0.47fc7 

0.439c 

0.4264 

0.4412 

0.1269 

0.1245 

0.1268 

0.4315 

0.4471 

0,4323 

0.5930 

0.397’ 

0. 3c5E 

0,405 7 

0.1237 

0.0953 

0,1211 

0.478c 

0.539. 

0.4752 

0.625c 

0.32<5 

{  7  Z  ' 

;  •  "  "  Z  Z 

0.1296 

0.1444 

0.1274 

0,545' 

0.5n7c 

0.5'" 

0.9419 

0.3252 

0.3349 

0.3236 

0.1567 

0.12c7 

0.1465 

0.5132 

0.5365 

0.52  =  7 

1 

p  • 

2.325c 

0.2524 

0.4174 

0.22-2 

r.  ?s.‘.4 

0.007' 

0.3353 

0.3972 

0.5866 

0.4356 

\; 

r. 

v  /[iElTA 

Uv'BUl 

UV/QM2 

uv/ Git* 

6AN*4 

1 

BANNA 

i 

r.' 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

uv 

Uv 

NV 

h. 

>9 

y>  _ 

ZCNE 

ZONE 

If 

II 

If 

II 

II 

II 

II 

M 

II 

II 

II 

II 

0. 127c 

-0.122c 

-0.116c 

-0.1237 

99.2160 

1.4007 

99.1680 

1.0005 

0 . 1 6 1 2 

-2.1255 

-0.1232 

-0.1254 

99.3720 

1.6008 

99.2750 

1.0015 

%  , 

0.2442 

-0.1173 

-0.1079 

-0.1164 

99.3160 

1.4007 

98.9860 

1.8009 

0.3023 

-0.1153 

-0.1244 

-0.1135 

98.9140 

3.4017 

98.5300 

6.6034 

0.3605 

-0.1134 

-0.1111 

-0,1129 

98.9160 

5.2027 

96.6880 

11.4060 

0.4767 

-0.09cE 

-0.079c 

-0.0959 

96.2470 

15.2080 

88.2220 

29.4150 

0.5930 

-0.0833 

-0.0432 

-0.0831 

92.3640 

24.4130 

75.8480 

46.4240 

0.8256 

-0.0459 

-0.0128 

-0.0497 

79,6440 

48.4250 

49.8160 

60.0310 

0.94JC 

-0.0401 

-0.0097 

-0.0423 

76.5830 

53.C270 

40.0690 

56.2290 

2.3256 

0.0101 

0.0000 

0,0111 

68.2740 

63.2320 

20.8480 

45.4230 
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Y/DELTA 

qll2/UEH2 

qll2/UEII2 

qlt2/UEII2 

UV/UEII2 

UV/UEU2 

UV/UEII2 

INTER- 

INTER- 

COMPOSITE 

TURBULENT 

TURBULENT 

COMPOSITE 

TURBULENT 

TURBULENT 

ZONE 

ZONE 

ZONE 

ZONE 

isssscssr 

:::ss::ss:iss::s:::ssss: 

xsssszssssss; 

0.1279 

0.011537 

0.011372 

0.011609 

-1.414E-03 

-1.349E-03 

-1.4J6E-03 

0.1860 

0.009976 

0.010099 

0,009977 

-1.252E-03 

-1.244E-03 

-1.252E-03 

0.2442 

0.009417 

0.009559 

0.009410 

-1.104E-03 

-1.032E-03 

-1.114E-03 

0.3023 

0.008465 

O.OOB267 

0.008547 

-9.758E-04 

-1.032E-03 

-9.697E-04 

0.3605 

0.007417 

0.006620 

0.007486 

-8.411E-04 

-7.353E-04 

-B.450E-04 

0 . 4  7  c  7 

0.006011 

0.004093 

0.006233 

-5.822E-04 

-3.260E-04 

-5.974E-04 

0.5930 

0.005376 

0.003961 

0.005714 

-4.478E-04 

-1.709E-04 

-4.746E-04 

0.8256 

0.004255 

0.003523 

0.004749 

-1.955E-04 

-4.509E-05 

-2.35’E-04 

0.941? 

0.003878 

0.003079 

0.004640 

-1.555E-04 

-2.981E-05 

-1 . 96 1 E "04 

2.325c 

0.003152 

0.001986 

0.003580 

3.192E-05 

0 . OOOE  *00 

3.985E-05 
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